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Abstract 


The requirement to settle or to position liquid fluid over the 
outlet end of a spacecraft propellant tank prior to main engine 
restart poses a microgravity fluid behavior problem. Resettlement 
or reorientation of liquid propellant can be accomplished by 
providing optimal acceleration to the spacecraft such that the 
propellant is reoriented over the tank outlet without any vapor 
entrainment, any excessive geysering, or any other undesirable fluid 
motion for the space fluid management under microgravity environment. 
The purpose of the present study is to investigate most efficient 
technique for propellant resettling through the minimization of 
propellant usage and weight penalties. Both full-scale and sub- 
scale liquid propellant tank of Space Transfer Vehicle have been used 
to simulate flow profiles for liquid hydrogen reorientation over the 
tank outlet. In sub-scale simulation, both constant and impulsive 
resettling acceleration have been used to simulate the liquid flow 
reorientation. Comparison between the constant reverse gravity 
acceleration and impulsive reverse gravity acceleration to be used 
for activation of propellant resettlement, it shows that impulsive 
reverse gravity thrust is superior to constant reverse gravity thrust 
for liquid reorientation in a reduced gravity environment. 
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(I) Statement of Problems 

Behavior of liquid propellant becomes fairly uncertain when the 
gravity environment reduces to the order of 10“ 6 g 0 level. The 
requirement to settle or to position liquid fuel over the outlet end of 
the spacecraft propellant tank prior to main engine restart poses a 
microgravity fluid behavior problem. Retromaneuvers of spacecraft, 
such as Orbital Maneuvering Vehicle (OMV) and Space Transfer Vehicle 
(STV) taking a flight from high earth orbit to low earth orbit, 
require settling or reorientation of propellant prior to main engine 
firing. Cryogenic liquid propellant is required to position over the 
tank outlet by using small auxiliary thrusters which provide a thrust 
parallel to the tank's major axis in the direction of flight. During 

the reorientation process, the liquid flows in an annular sheet along 

\ 

the tank wall, with the gas or ullage bubble "rising" centrally into 
the liquid. This motion would also clear the tank vent of liquid so 
that venting of vapor is possible. 

An efficient propellant settling technique should minimize 
propellant usage and weight penalties. This can be accomplished by 
providing optimal acceleration to the spacecraft such that the 
propellant is reoriented over the tank outlet without any vapor 
entrainment, any excessive geysering, or any other undesirable fluid 
motion. 

Time-dependent dynamical behaviors of fluids in a microgravity 
environment have been investigated by numerically computing full set 
of Navier-Stokes equation 1 2 “ 1 5 1 . 
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(II) Problems Attempted 


Liquid propellant (hydrogen) tank is 168 inch in diameter, has a 
square root of 2 domes and a 48 inch long barrel section. Figure 1 
shows the size and dimension of hydrogen tank. A settling 
acceleration, ranging from 1.5 x 10“ 3 to 3.0 x 10 3 g 0 , has been 
assumed to be produced by idle-mode thrust from the main engine. 
Settling of propellant prior to main engine start is defined to be 
required at two tanks with fill percentages of 30% and 70%. 

To provide conservative settling conditions, the liquid is 
assumed to be initially oriented at the top of the tank. If the 
spacecraft has been coasting for a long time, aligned with its 
direction of motion, the most significant force, drag would be axial. 

Axisymmetric representation of the tank, with a cylindrical 
coordinate system is assumed. The settling acceleration is assumed 
to be parallel to by the axial direction, and the liquid is initially 
oriented symmetrically at the top of the tank. 

Numerical simulations have been carried out for both full-scale 
and sub-scale of liquid hydrogen propellant tanks for the Space 

Transfer Vehicle (STV) . 

(Ill) Full-Scale Simulation for Liquid 
Reorientation and Resettlement 

The following four cases have been accomplished for the cases of 
full-scale STV propellant tank liquid fluid reorientation: (A) 

Liquid filled level of 30% and settling acceleration of 1.5 x 10 g 0 ; 
(B) Liquid filled level of 30% and settling acceleration of 3.0 x 10 
g 0 ; (C) Liquid filled level of 70% and settling acceleration of 1.5 x 
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Figure 1 Size and Diameter of Spa 


10" 3 g 0 ; and (D) Liquid filled level of 70% and settling acceleration 
of 3.0 x 10 -3 g 0 . , 

Figures 2 to 11 show the time evolution of liquid hydrogen 
reorientation for liquid filled level of 30% and settling 
acceleration of 1.5 x 10" 3 g 0 . Figure 12 to 21 illustrate the time 
evolution of liquid hydrogen reorientation for liquid filled level of 
30% and resettling accleration of 3.0 x 10 3 g 0 . Figure 22 to 31 show 
the time evolution of liquid hydrogen reorientation of liquid 
hydrogen for liquid filled level of 70% and resettling acceleration of 
1.5 x 10' 3 g 0 . Figure 32 to 41 illustrate the time evolution of liquid 
hydrogen for liquid filled level of 70% and resettling of 3.0 x 10 

9o- 


5 


(A) Full-Scale Simulation 
Liquid Filled: 30% 

Settling Acceleration: 1.5 x 10 3 g 0 
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(B) Full-Scale Simulation 
Liquid Filled: 30% 

Settling Acceleration: 3.0 x 10“ 3 g 0 
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Figure 13 
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Figure 18 
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Figure 20 



CN 

CD 


CD 

Ll_ 



I ' i : I r i ' I 

‘OOfr O’OOC Q‘002 0*001 0*0 


(UI0 ) IH9I3H 


27 


RADIUS (CM) 



(C) Full-Scale Simulation 
Liquid Filled: 70% 

Settling Acceleration: 1.5 x 10" 3 g 0 
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Figure 22 
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Liquid Hydrogen and Vapor 
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(D) Full-Scale Simulation 
Liquid Filled: 70% 

Settlinng Accleration: 3.0 x 10" 3 g 0 
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Figure 36 


Liquid hydrogen and Vapor 
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Figure 39 
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Figure 40 
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(IV) Sub-Scale Simulation for Liquid 
Reorientation and Resettlement 

In the sub-scale STV propellant tank liquid hydrogen 
reorientation we have accomplished the following four cases of 
numerical simulation: (E) Cryogenic liquid hydrogen reoreintation 
activated by constant settling acceleration of geyser initiation; (F) 
Cryogenic liquid hydrogen reorientation activated by low frequency 
impulsive settling acceleration of geyser initiation; (G) Cryogenic 
liquid hydrogen reorientation activated by medium frequency 
impulsive settling acceleration of geyser initiation; and (H) 
Cryogenic liquid hydrogen reorientation activated by high frequency 
impulsive settling acceleration of geyser initiation. 
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(E) Cryogenic Liquid Hydrogen Reorientation 
Activated by Constant Settling Acceleration 
of Geyser Initiation 
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Abstract 


A key objective of the cryogenic fluid management of the 
spacecraft propulsion system is to develop the technology necessary 
for acquisition or positioning of liquid outflow or vapor venting. 
In this paper, numerical simulation of positive liquid acquisition is 
attempted by introducing reverse gravity acceleration, resulting 
from the propulsive thrust of auxiliary engines, which exceeds 
critical value for the initiation of geyser. Based on the computer 
simulation of flow fields during the course of fluid reorientation, 
six dimensionless parameters resulted in this study. It shows that 
these parameters hold near constant values through the entire ranges 
of liquid filled levels, from 3 0% to 80%, during the course of fluid 
reorientation . 
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Nomenclature 


* 

a g = geyser initiation acceleration (cm/s ) 

a = scale flow acceleration associated with maximum velocity 
(cm/s 2 ), defined by Eq (5) 

D = diameter of propellant tank (cm) 
f = frequency of impulsive thrust (H z ) 
g i = geyser initiation gravity-level (g 0 ) 

2 

g 0 = normal Earth gravitational acceleration = 9.81 m/s 
h = average liquid height (cm) 

h ra = maximum liquid height (cm) 

L = height of propellant tank (cm) 

L = scale length of maximum liquid height (cm/s) , defined by 
m 

Eq (4) 

STV = Space Transfer Vehicle 
t f = average free fall time (s) 

t = time for observing maximum flow velocity (s) 
t R = liquid reaching tank bottom time (s) 

u — free fall velocity from maximum liquid height (cm/s) , defined by 

V f HI 

Eq (3) 

V f = average free fall velocity (cm/s), defined by Eq (2) 

V = maximum flow velocity (cm/s) 
m 
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I. Introduction 


In spacecraft design, the requirements fora settled propellant 
are different for tank pressurization, engine restart, venting, or 
propellant transfer. Prepressurization requires that heat and mass 
transfer effects be minimized; otherwise, a process of chill down of 
tank, venting of noncondensing gases, etc. , may have to carry out for 
the cryogenic system. For engine restart, it is necessary to have the 
liquid settle with no bubbles near the tank outlet so that the initial 
flow of propellant will not carry vapor to the pump or engine. The 
slosh wave amplitude should be relatively low to keep the center of 
mass shifts within an acceptable range and wave motion low enough to 
avoid pressure collapse caused by interface agitation . 

Bubble and globule formation as a result of liquid impact with 
the aft end of the tank could lead to propellant loss for the 
spacecraft during venting. Globules could be entrained in the vented 
ullage gas or bubbles rising through the liquid and expanding because 
of the decreasing tank pressure could cause a spray of globules to be 
vented. Liquid level rise, vent liquid loss, fluid freezing, and 
vehicle dynamics are all affected by the microgravity levels. 

A key objective of the cryogenic fluid management of spacecraft 
propulsion system, such as a Space Transfer Vehicle 1 (STV) , is to 
develop the technology necessary for acquisition or positioning of 
liquid and vapor within a tank in reduced gravity to enable liquid 
outflow or vapor venting. Liquid acquisition techniques can be 
divided into two general categories: (1) Active liquid acquisition 
by the creation of a positive acceleration environment resulting from 
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the propulsive thrust of small auxiliary engines, and (2) Passive 
liquid acquisition utilizing the liquid capillary forces provided by 
using solid baffles of'liquid traps made of fine mesh screen material. 
In this study, active liquid acquisition is aimed for numerically 
simulating the resettlement of cryogenic liquid hydrogen. Liquid 
hydrogen, which, in general, poses more severe technical challenges 
than liquid oxygen, is used as the test bed working fluid in this 

study . 

Recently Lesie 2 was able to measure and to numerically compute 

the bubble shapes at various ratios of centrifugal force to surface 

, ^ 

tension force in the microgravity environment. Hung and Leslie 
extended Leslie's work 2 to rotating free surfaces influenced by 
gravity with higher rotating speeds when the bubble intersects with 
both the top and bottom walls of the cylinder. Hung et al. further 
extended the work to include rotating speeds which resulted with 
bubbles intersecting and/or without intersecting the top, bottom and 
side walls of the cylinder. 

An analysis of time-dependent dynamical behavior of surface 
tension on partially-filled rotating fluids in both low gravity and 
microgravity environments was carried out by numerically solving the 
Navier-Stokes equations subjected to the initial and the boundary 
conditions 4 ' 6 . The initial condition for the bubble profiles was 
adopted from steady-state formulations developed by Hung and Leslie 3 , 
and Hung et al. 6 for rotating cylinder tank; and by Hung et al. for 
the dewar-shaped container to be used in the Gravity Probe-B 
Spacecraft 8 . Some of the steady-state formulations of bubble 
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shapes, in particular for bubbles intersecting at the top wall of the 
cylinder, were compared with the experiment carried out by Leslie in 
a free-fallinng aircraft (KC-135) with an excellent agreement. 

An efficient propellant settling technigue should minimize 
propellant usage and weight penalties. This can be accomplished by 
providing optimal acceleration to the spacecraft such that the 
propellant is reoriented over the tank outlet without any vapor 
entrainment, any excessive geyser ing, or any other undesirable fluid 

motion. 

Production of geyser during the propellant reorientation is not 
a desirable motion for the space fluid management. It is because 
geyser is always accompanied by the vapor entrainment and globule 
formation. Geyser is observed at reverse gravity acceleration 
greater than certain critical values of acceleration during the 
course of liquid reorientation. In other words, geyser will not be 
observed at very low reverse gravity level, and it will be detected 
when the reverse gravity level is greater than the certain critical 
level. In this paper, numerical simulation of positive liquid 
acquisition is attemted by introducing reverse gravity acceleration, 
resulting from the propulsive thrust of small auxiliary engines which 
exceeds the critical value for geyser initiation. The reverse 
gravity acceleration is starting with a small value and increases 
gradually till the initiation of geyser is detected in the computer 
simulation for the liquid reorientation of propellant tank with 
various liquid— f illed levels. 

In this study, time-dependent computations have been carried out 
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to investigate the dynamical behavior of fluid reorientation or 
resettling of propellant prior to main engine firing for spacecraft 
restart at the net reverse gravity acceleration which is great enough 
to initiate geyser during the liquid reorientation. The computation 
extends to the study of the characteristics of fluid resettlement due 
to the reversal of lowest artificial gravity fields which is great 
enough to initiate geyser with and without various frequencies of 
impulsive acceleration. This frequency of impulsive acceleration is 
generally termed "gravity j itters" . Gravity j itters are produced by 
spacecraft attitude motion, machinery (turbine, pump, engine) 
vibrations, thruster firing, thruster shutdown, etc . 9 Positioning 
of liquid propellant over the tank outlet can be carried out by using 
small auxiliary thrusters which provide a thrust parallel to the 

tank's major axis in the direction of flight. 

H. Numerical simulation of Liquid Hydorgen Reorientatio n at 
the Reverse Gravity Acceleration of Geyser Initiatio n 

Propellant tank of liquid-filled levels of 30, 50, 65, 70 and 80% 
are considered in this study. Time-dependent axial symmetry 
mathematical formulation are adopted. Detailed description of 
mathematical formulation, initial and boundary conditions suitable 
for the analysis of cryogenic fluid management under microgravity 
environment are given in our earlier studies . 4 , 6 , 10-12 The initial 
profiles of liquid-vapor interface are determined from computations 
based on algorithms developed for the steady state formultion of 

3 _ 7 

microgravity fluid management. 

A staggered grid for the velocity components is used in this 
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computer program. The method was first developed by Harlow and 
Welch 13 for their MAC (marker and cell) method of studying fluid flows 
along free surface. The finite difference method employed in this 
numerical study was the "Hybrid Scheme" developed by Spalding . 1 4 The 
formulation for this method is valid for any arbitrary interface 
location between the grid points is not limited to middle point 
interfaces. 15 An algorithm for semi-implicit method 16 was used as 
the procedure for modeling the flow field. The time step is 
determined automatically based on the size of the grid points and the 
velocity of flow fields. Detailed description of the computational 
algorithm applicable to microgravity fluid management are 
illustrated in our earlier studies. 4 ' 6 ' 1 °“ 1 2 Fig 1(A) shows the 
distribution of grid points in the radial-axial plane of cylindrical 

coordinates. 

For the purpose of facilitating easy comparison between 
computational resuslts and experimental measurement, a model of 0.01 
size prototype is adopted in the computer simulation. Fig 1(B) shows 
a model size for computation. The size of prototype is height, L - 
166.8" (423.672 cm), and diameter, D= 168" (426.72 cm). Model size 
is L = 4.23672 cm and D = 4.2672 cm, as shown in Fig 1(B) . If the 
spacecraft had been coasting for a long time, aligned with its 
direction of motion, the most signigicant force, drag, would be axial 
and with acceleration of 10“ 4 g 0 along upward direction. The hydrogen 
vapor is, thus, originally positioned at the bottom of the tank. The 
requirement to settle or to position liquid fuel over the outlet end of 
the spacecraft propellant tank prior to main engine restart poses a 
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Retromaneuvers of 


microgravity fluid behavior problem. 10 
spacecraft, such as STV, require settling or reorientation of the 
propellant to main engine firing. 1 ' 10 Cryogenic liquid propellant 
shall be positioned over the tank outlet by using auxiliary thrusters 
(or idle-mode thrusters from the main engine) which provide a thrust 
parallel to the tank's major axis in the direction of flight. In the 
present study of computer simulation, a small value of reverse gravity 
acceleration (downward direction) is provided by the propulsive 
thrust of small auxiliary engine to initiate the reorientation of 
liquid propellant. This small value of reverse gravity acceleration 
of propulsive thrust increases gradually till reaching the critical 
value on which initiation of geyser is detected during the time period 
of fluid resettlement. We term this reverse gravity acceleration of 
propulsive thrust, which is capable to initiate geyser, as "geyser 
initiation gravity -level" . This geyser initiation gravity level has 
been investigated through the method of trial and error for the 
various liquid-filled levels of propellant tank as a base to simulate 
following cases of reduced gravity fluid behaviors during the 
reorientation: (A) constant reverse gravity acceleration, (B) 

impulsive reverse gravity acceleration with a frequency of 0.1 H z , (C) 
impulsive reverse gravity acceleration with a frequency of 1. 0 H z , and 
(D) impulsive reverse gravity acceleration with a frequency of 10 
H Cases (B) to (D) for the impulsive reverse gravity acceleration 
with various frequencies will be discussed in the subsequent papers. 
Liquid filled level of 30, 50, 65, 70 and 80% have been considered in 
this study. Cryogenic liquid hydrogen at temperature of 20K is 
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considered. Hydrogen density of 0.071 g/cm 3 ; surface tension 
coefficient at the interface between liquid hydrogen and hydrogen 
vapor of 1.9 dyne/cm; hydrogen viscosity coefficient of 1.873 x 10" 3 
cm 2 /s; and contact angle of 0.5° are used in the computer simulation. 

Reorientation of cryogenic liquid hydrogen activated by geyser 
initiation reverse gravity acceleration produced by propulsive 
thrust has been investigated for various liquid filled levels of 
propellant tank. It is found that these geyser initiation gravity 
levels are 5.5 x 10“ 2 , 6.52 x 10' 2 , 6.6 x 10‘ 2 , 6.7 x 10“ 2 and 8.2 x 
10“ 2 g 0 for liquid filled levels of 30, 50, 65, 70, and 80%, 
respectively. To illustrate some examples Figs 2 and 3 show the 
selected sequences of time evolution of fluid reorientation for 
cryogenic hydrogen with liquid filled levels of 30, and 80%, 
respectively. Each figure contains four sub— figures. Subfigure 
(A) is initial profile of liquid-vapor interface at the moment of the 
starting of fluid reorientation at time t = 0; subfigure (B) , the flow 
profile during the course of fluid reorientation before the initation 
of geysering motion; subfigure (C) , the flow profile with geysering 
motion; and subfigure (D) , the flow profile after the ending of 
geysering motion. 

Figs 2 and 3 illustrate following flow behaviors: (1) The 
liquid starts to flow in an annular sheet along the solid wall of tank 
and gradually pushes the vapor toward the central portion of the lower 
dome of tank as the net acceleration, reversing the direction of 
gravity field, which is applied toward the downward direction of 
tank's major axis, by using small auxiliary thrusters; (2) As the 
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downward fluid annular sheet along the tank wall reaches the central 
bottom dome side of the tank, a geysering flow is observed? and (3) The 
vapor is thus pushed upward centrally into the liquid and the 
geysering disappears . 

Table 1 shows the characteristics of cryogenic liquid hydrogen 
resettlement activated by reverse gravity acceleration at geyser 
initiation gravity level. Average liquid height h, and maximum 
liquid height h m are shown in Fig 1(B). Average free fall velocity 
V f , average free fall time t f , and free fall velocity from maximum 
liquid height V fm are computed from the following equations: 


Vf = 

(2g t h) 1/2 

(1) 


( 2h l 

( 1 / 2 

(2) 


lg J 



Vfm = 

(2g 

•h ) 1/2 
1 1A m > 

(3) 


where q i denotes geyser initiation reverse gravity acceleration. 
Values of maximum flow velocity V m , time for observing maximum flow 
velocity t m , and time for reorienting liquid flowing down and reaching 
the bottom of propellant tank t R , are obtained from the numerical 
computation of flow field. Scale length of maximum flow velocity L m , 
and scale flow acceleration associated with maximum velocity a m , are 
computed from the following parameters: 

L_ = V m t_ (4) 

XJ Tn m m 
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(5) 



Following dimensionless parameters are introduced: V m /V f/ t R /t f , 
t m /t f , a m /a g , L m /h and V m /V fra where a g stands geyser initiation 
acceleration (cm/s 2 ) for corresponding geyser initiation gravity 
level g i . 

Figs 4 to 6 show the variations of dimensionless parameters in 
terms of liquid filled levels. Denominators of these six 
dimensionless parameters are either predetermined from the geometry 
of liquid fill levels or can be deduced from the corresponding 
calculations associated with the geyser initiation gravity levels. 
Characteristics of these near constant values dimensionless 
parameters can provide a good understanding of the physics of 
microgravity fluid behaviors, in particular the active category of 
liquid acquisition or positioning, and also the design criteria of on- 
orbit spacecraft propulsion system at the critical value of reverse 
gravity acceleration of propulsive thrust which is capable to 
inititiate geyser. 

Fig 4 (A) shows the ratio of maximum flow velocity to average free 
fall flow velocity V m /V f and its associated parameters of V m and V f in 
terms of liquid filled levels. It shows that the ratio of V m /V f 
varies in the range of 4.3 to 4 . 9 in the entire liquid filled levels 
while V m and V f vary from 62.0 to 73.6 cm/ s (decreasing with increasing 
liquid filled levels) and from 12.5 to 17.0 cm/s (also decreasing with 
increasing liquid filled levels), respectively. As V f can be 
predetermined from geyser initiation gravity level and average liquid 
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height, shown in Eq (1), one can make an approximate prediction of 
maximum flow velocity during the liquid reorientation for the various 
liquid filled levels. 

Fig 4(B) shows the ratio of liquid reaching bottom time to 
average free fall time t R /t f and its associated parameters of t R and t f 
in terms of liquid filled levels. It shows that the ratio of t R /t f 
varies in the range of 1.21 to 1.3 0 in the entire liquid filled levels 
while t R and t f vary from 0.20 to 0.40 s (decreasing with increasing 
liquid filled levels) and from 0.15 to 0.32 s (also decreasing with 
increasing liquid filled levels) , respectively. As t f can be 
predetermined from geyser initiation gravity level and average liquid 
height, shown in Eq (2) , one can predict the time reorienting liquid 
fluid flowing down from the original position and reaching the bottom 
of propellant tank for the various liquid filled levels at the reverse 
gravity acceleration capable for the initiation of geyser. 

Fig 5(A) shows the ratio of time for observing maximum flow 
velocity to average free fall time t m /t f and its associated parameters 
of t ra and t f in terms of liquid filled levels. It shows that the ratio 
of t TO /t f varies in the range of 1.2 to 1.3 in the entire liquid filled 
levels while t m and t f vary from 0.20 to 0.42 s (decreasing with 
increasing liquid filled levels) and from 0.15 to 0.32 s (also 
decreasing with increasing liquid filled levels) , resectively. As 
we indicated in 4, t f can be predetermined, one can predict the time 
for observing maximum flow velocity for various liquid filled levels 
at the reverse gravity acceleration capable for the initiation of 

geyser. 
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Fig 5(B) shows the ratio of scale flow acceleration associated 
with maximum velocity to geyser initiation acceleration for 
corresponding geyser initiation gravity level a^/ag and its 
associated parameters of a m and a g in terms of liquid filled levels. 
It shows that the ratio of a m /a g varies in the range of 3 . 3 to 3 . 9 in the 
entire liquid filled levels while a m and a g vary from 175 to 310 cm/s 2 
(increasing with increasing liquid filled levels) and from 53.9 to 
80.4 cm/s 2 (also increasing with increasing liquid filled levels) , 
respectively. As a g can be predetermined from geyser initiation 
gravity level, one can make an approximate prediction of scale flow 
acceleration associated with maximum velocity, which is defined in Eq 
(5) , at the reverse gravity acceleration capable for the initiation of 

geyser. 

Fig 6 (A) shows the ratio of scale length of maximum flow velocity 
to average liquid height L ro /h and its associated parameters of L m and h 
in terms of liquid filled levels. It shows that the ratio of L m /h 
varies in the range of 11.6 to 13.0 in the entire liquid filled levels 
while L m and h vary from 12.4 to 30.9 cm (decreasing with increasing 
liquid filled levels) and from 0.93 to 2.67 cm (also decreasing with 
increasing liquid filled levels), respectively. As h can be 
predetermined from the geometry of liquid filled levels, one can make 
an approximate prediction of scale length of maximum flow velocity, 
which is defined in Eq (4), at the reverse gravity acceleration 

capable for the initiation of geyser. 

Fig 6(B) shows the ratio of maximum flow velocity to free fall 
velocity from maximum liquid height V m /V fm and its associated 
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parameters of V m and V fm in terms of liquid filled levels. It shows 
that the ratio of V m /V f m varies in the range of 3 . 8 to 4 . 0 in the entire 
liquid filled levels 'while V m and V fm vary from 62.0 to 73.6 cm/s 
(decreasing with increasing liquid filled levels) and from 16.4 to 
19.2 cm/s (also decreasing with increasing liquid filled levels), 
respectively. As V fra can be predetermined from geyser initiation 
gravity level and maximum liquid height, shown in Eq (3) , one can make 
an approximate prediction of maximum flow velocity at the reverse 
gravity acceleration capable for the initiation of geyser. 

Six dimensionless parameters presented in this study show that 
the parameters hold near constant values through the entire ranges of 
liquid filled levels during the course of reorientation of liquid 
hydrogen activated by the reverse gravity acceleration which is great 
enough to initiate geysering flow. As the denominators of these six 
dimensionless parameters are either predetermined from the geometry 
of liquid filled levels or can be deduced from the corresponding 
calculations associated with the geyser initiation gravity levels, 
one can predict the flow parameters from these relations. 

Ill . Discussion and Conclusions 

The efficient management of subcritical cryogenic propellants 
is one of the key technology drivers for the on-orbit spacecraft. 
Cryogenic liquids are essential for the spacecraft as reactants, 
coolants, and propellants. The requirement to settle or to position 
liquid fuel over the outlet end of the spacecraft propellant tank 
prior to main engine restart poses a microgravity fluid behavior 
problem. Retromaneuvers of spacecraft require settling or 
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ion of the propellant prior to main ®ngine firing* 
Cryogenic liquid propellant is positioned over the tank outlet by 
using small auxiliary thrusters (or idle-mode thrusters from the main 
engine) which provide a thrust parallel to the tank's major axis in the 
direction of flight. 

The results of the study of fluid reorientation have to be 
evaluated in terms of how well they can be managed efficiently. An 
efficient propellant settling technique should minimize propellant 
usage and weight penalties through the operation of small thrusters 
(or idle— mode thrusters from the main engine) . This can be 
accomplished by providing optimal acceleration to the spacecraft such 
that the propellant is reoriented over the tank outlet without any 
vapor entrainment, any excessive geysering, or any other undesirable 
fluid motion. 

Production of geyser during the propellant reorientation is not 
a desirable motion for the space fluid management under microgravity 
environment. In this paper, numerical simulation of positive liquid 
acquisition is attempted by introducing reverse gravity 
acceleration, resulting from the propulsive thrust of auxiliary 
engines, which exceeds critical value for the intiation of geyser. 

Flow profile simulations, in particular the time evolution of 
liquid-vapor interface during the course of fluid reorientation, for 
the selected various liquid filled levels are shown in Figs 2 and 3. 
Table 1 and Figs 4 to 6 show the characteristics of dimensionless flow 
parameters and their associated flow fields during the course of fluid 
reorientation with geyser initiation. Computer simulation of flow 
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fields disclose the following results: (1) Geyser initiation 
gravity level is high (low) for high (low) liquid filled level? (2) 
Average liquid height is high (low) for low (high) liquid filled 
level; (3) Average free fall flow velocity is high (low) for low 
(high) liquid filled level; (4) Maximum flow velocity is high (low) 
for low (high) liquid filled level; (5) Average free fall time is 
longer (shorter) for low (high) liquid filled level; (6) Time for 
liquid falling from original position to liquid reaching bottom time 
is longer (shorter) for low (high) liquid filled level ; (7) Time for 
observing maximum flow velocity is longer (shorter) for low (high) 
liquid filled level; (8) Scale length of maximum flow velocity is 
longer (shorter) for low (high) liquid filled level; (9) Scale flow 
acceleration associated with maximum velocity is low (high) for low 
(high) liquid filled level; (10) Maximum liquid height is high (low) 
for low (high) liquid filled level; and (11) Free fall velocity from 
maximum liquid height is high (low) for low (high) liquid filled 
level . 

Based on the computer simulation of flow fields during the course 
of fluid reorientation, six dimensionless parameters are presented in 
this study. It is shown, in Figs 4 to 6, that these parameters hold 
near constant values through the entire ranges of liquid filled levels 
during the course of fluid reorientation activated by the reverse 
gravity acceleration great enough to initiate geyser. As the 
denominators of these dimensionless parameters are either 
predetermined from the geometry of liquid filled levels or can be 
deduced from the corresponding calculations associated with the 
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geyser intiation gravity levels, one can predict the values of these 
flow parameters from the relationship shown in Eqs (1) to (5) . 

Any fluid capable of motion relative to the spacecraft will be 
subject to an acceleration relative to the mass center of the 
spacecraft that arises from the gravity gradient of the Earth 17 ' 18 . 
In addition to the Earth's gravitational force, the interaction 
between the particle mass of fluids and the spacecraft mass due to 
gravity gradient accelerations 1 7 have also been taken into 
consideration in this microgravity fluid management study. 

To conclude, we have demonstrated that, the computer algorithm 
presented, can be used to simulate fluid behavior in a microgravity 
environment, in particular the development of technology necessary 
for acquisition or positioning of liquid and vapor within a tank to 
enable liquid outflow or vapor venting through active liquid 
acquisition by the creation of a positive acceleration environment 
resulting from propulsive thrust. Better understanding of the full 
pictures of flow fields during the course of fluid reorientation can 
provide the proper design techniques for handling and managing the 
cryogenic liquid propellants to be used in on-orbit spacecraft 
propulsion. 
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Table 1 CHARACTERISTICS OF CRYOGENIC HYDROGEN REORIENTATION 
(Frequency of Impulsive Acceleration » 0.0 Hz) 



Liquid Reaching Bottom Time, t R (s) 0.40 0.31 0.27 0.25 0.20 
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Figure Captions 


Fig 1. 


Fig 2 . 


Fig 3 . 


Fig 4 . 


Fig 5. 


Fig 6 . 


(A) Distribution of grid points in the radial— axial plane 
of cylindrical coordinate for propellant tank, and (B) 
Model size propellant tank adopted for numerical 
simulation with geometrical description. 

Selected sequences of time evolution of fluid 
reorientation with liquid filled level of 3 0% , (A) initial 
profile, (B) flow prof ile before the initiation of geyser , 
(C) flow profile with geyser, and (D) flow profile after 
the ending of geyser. 

Selected sequences of time evolution of fluid 
reorientation with liquid filled levels of 80-3, (A) 

Initial profile, (B) flow profile before the initiation of 
geyser, (C) flow profile with geyser, and (D) flow profile 
after the ending of geyser. 

(A) Ratio of V m /V f and its associated parameters in terms 
of liquid filled levels, (B) Ratio of t R /t f and its 
associated parameters in terms of liquid filled levels. 
(A) Ratio of t m /t f and its associated parameters in terms 
of liquid filled levels, (B) Ratio of a m /a g and its 
associated parameters in terms of liquid filled levels. 
(A) Ratio of L m /h and its associated parameters in terms of 
liquid filled levels, (B) Ratio of V m /V fra and its 
associated parameters in terms of liquid filled levels. 
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ABSTRACT 


The requirement to settle or to position liquid fluid over the 
outlet end of spacecraft propellant tank prior to main engine restart 
poses a microgravity fluid behavior problem. Resettlement or 
reorintation of liquid propellant can be accomplished by providing 
optimal acceleration to the spacecraft such that the propellant is 
reorientaed over the tank oulet without any vapor entrainment, any 
excessive geysering, or any other undesirable fluid motion for the 
space fluid management under microgravity environment. The purpose 
of present study is to investigate most efficient technique for 
propellant resettling through the minimization of propellant usage 
and weight penalties. Comparison between the constant reverse 
gravtiy acceleration and impulsive reverse gravity acceleration to be 
used for the activation of propellant resettlement, it shows that 
impulsive reverse gravity thrust is superior to constant reverse 
gravity thrust for liquid reorientation in a reduced gravity 


environment. 


I . Introduction 


In spacecraft design, the requirements fora settled propellant 
are different for tank pressurization, engine restart, venting, or 
propellant transfer. Prepressurization requires that heat and mass 
transfer effects be minimized; otherwise, a process of chill down of 
tank, venting of noncondensing gases, etc. , may have to carry out for 
the cryogenic system. For engine restart, it is necessary to have the 
liquid settle with no bubbles near the tank outlet so that the initial 
flow of propellant will not carry vapor to the pump or engine. The 
slosh wave amplitude should be relatively low to keep the center of 
mass shifts within an acceptable range and wave motion low enough to 
avoid pressure collapse caused by interface agitation . For venting, 
it is probably necessary that virtually all bubbles be displaced from 
the bulk liquid so that a two-phase mixture is not vented. Propellant 
transfer requires that the liquid be completely settled with 
virtually no bubbles. Outflow of a liquid near the tank outlet can 
result in the prematrue ingestion of gas while a signigicant amount of 
liquid is still in the tank under microgravity environment. This 
phenomenon is termed "suction dip". Slosh wave motion must be 
minimal because the combination of "suction dip" and sloshing could 
cause gas pull-through to occur more readily in microgravity than if 
the surface were essentially quiescent. 

During the prepressurization of a cryogenic propellant in 
microgravity, significant heat and mass transfer will occur if the 
liquid interface is disturbed. Interface disturbances may result 
from (a) impingement of the gas on the liquid surface at a mass flow 
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rate sufficient to cause Kelvin-Helmholtz instability, (b) globule 
formation from breaking waves caused by wave motion over baffles or 
internal hardware., (c) globule and surface froth formation resulting 
from movement of bubbles through the liquid to the surface, and (d) 
surface froth formation because of gas impingement. 

Bubble and globule formation as a result of liquid impact with 
the aft end of the tank could lead to propellant loss for the 
spacecraft during venting. Globules could be entrained in the vented 
ullage gas or bubbles rising through the liquid and expanding because 
of the decreasing tank pressure could cause a spray of globules to be 
vented. Liquid level rise, vent liquid loss, fluid freezing, and 
vehicle dynamics are all affected by the microgravity levels. 

A key objective of the cryogenic fluid management of spacecraft 
propulsion system, such as a Space Transfer Vehicle 1 (STV) , is to 
develop the technology necessary for acquisition or positioning of 
liquid and vapor within a tank in reduced gravity to enable liquid 
outflow or vapor venting. Liquid acquisition techniques can be 
divided into two general categories: (1) Active liquid acquisition 
by the creation of a positive acceleration environment resulting from 
the propulsive thrust of small auxiliary engines, and (2) Passive 
liquid acquisition utilizing the liquid capillary forces provided by 
using solid baffles of liquid traps made of fine mesh screen material . 
In this series of study, active liquid acquisition is aimed for 
numerically simulating the resettlement of cryogenic liquid 
hydrogen. Liquid hydrogen, which, in general, poses more severe 
technical challenges than liquid oxygen, is used as the test bed 
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working fluid in this study . 

Recently Lesie 2 was able to measure and to numerically compute 
the bubble shapes at various ratios of centrifugal force to surface 
tension force in 2 , 4 and 6.3 cm deep cylinders in the microgravity 
environment. The results showed excellent agreement between model 
computation and measurements. Hung and Leslie 3 extended Leslie's 
work 2 to rotating free surfaces influenced by gravity with higher 
rotating speeds when the bubble intersects with both the top and 
bottom walls of the cylinder. Hung et al. 4 ' 5 further extended the 
work to include rotating speeds which resulted with bubbles 
intersecting and/or without intersecting the top, bottom and side 
walls of the cylinder. 

An analysis of time-dependent dynamical behavior of surface 
tension on partially-filled rotating fluids in both low gravity and 
microgravity environments was carried out by numerically solving the 
Navier-Stokes equations subjected to the initial and the boundary 
conditions 4 ' 6 . At the interface between the liquid and the gaseous 
fluids, both the kinematic surface boundary condition, and the 
interface stress conditions for components tangential and normal to 
the interface, were applied. The initial condition for the bubble 
profiles was adopted from steady-state formulations developed by Hung 
and Leslie 3 , and Hung et al . 6 for rotating cylinder tank; and by Hung 
et al. 5 ' 7 for the dewar-shaped container to be used in the Gravity 
Probe-B Spacecraft 8 . Some of the steady-state formulations of 
bubble shapes, in particular for bubbles intersecting at the top wall 
of the cylinder, were compared with the experiment carried out by 
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Leslie 2 in a free-fallinng aircraft (KC-135) . Comparisons of time- 
dependent results between numerical computations and experiments 
were unavailable. . This was because the calibration of the recordings 
of time-dependent gravity variations in a KC-135 aircraft during the 
short time periods of microgravity environment is very difficult. 
There was also an unavailability of accelerometer data for measuring 
the actual levels of microgravity during the experiment. 

An efficient propellant settling technique should minimize 
propellant usage and weight penalties. This can be accomplished by 
providing optimal acceleration to the spacecraft such that the 
propellant is reoriented over the tank outlet without any vapor 
entrainment, any excessive geyser ing, or any other undesirable fluid 
motion. 

Production of geyser during the propellant reorientation is not 
a desirable motion for the space fluid management. It is because 
geyser is always accompanied by the vapor entrainment and globule 
formation. Geyser is observed at reverse gravity thrust greater than 
certain critical values of acceleration during the course of liquid 
reorientation. In other words, geyser will not be observed at very 
low reverse gravity level, and it will be detected when the reverse 
gravity level is greater than the certain critical value. In this 
series of study, numerical simulation of positive liquid acquisition 
is attemted by introducing reverse gravity acceleration, resulting 
from the propulsive thrust of small auxiliary engines which exceeds 
the critical value for geyser initiation. The reverse gravity 
acceleration is starting with a small value and increases gradually 
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■j. in the initiation of geyser is detected in the computer simulation 
for the liquid reorientation of propellant tank with various liquid- 
filled levels. 

In this series of studies, time-dependent computations have been 
carried out to investigate the dynamical behavior of fluid 
reorientation or resettling of propellant prior to main engine firing 
for spacecraft restart at the net reverse gravity acceleration which 
is great enough to initiate geyser during the liquid reorientation. 
First paper of present study (Paper I ) 9 investigates the 
characteristics of fluid resettlement due to the reversal of lowest 
constant reverse gravity acceleration which is great enough to 
initiate geyser without any impulsive acceleration. The frequency 
of impulsive acceleration is generally termed "gravity jitters". 
Gravity jitters are produced by spacecraft attitude motion, machinery 
(turbine, pump, engine) vibrations, thruster firing, thruster 
shutdown, impulsive engine acceleration, etc . 10 Positioning of 
liquid propellant over the tank outlet can be carried out by using 
small auxiliary thrusters which provide a thrust parallel to the 
tank's major axis in the direction of flight. 

Computer simulation of flow field based on Paper I 9 during the 
course of fluid reorientation induced by constant reverse gravity 
acceleration show that six dimensionless parameters are abtained in 
the study. These parameters hold near constant values through the 
entire ranges of liquid filled levels during the course of fluid 
reorientation activated by the reverse gravity acceleration great 
enough to initiate geyser. As the denominators of these 
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dimensionless parameters are either predetermined from the geometry 
of liquid filled levels or can be deduced from the corresponding 
calculations associated with the geyser initiation gravity levels. 
One can predict the values of these flow parameters. These 
predictable parameters include maximum flow velocity V m , time for 
observing maximum flow velocity t m , time for reorienting liquid 
flowing down and reaching the bottom of propellant tank t R , scale 
length of maximum flow velocity L m , and scale flow acceleration 

associated with maximum velocity a B . 

Instead of applying constant reverse gravity acceleration as we 
described in Paper I 9 , this paper adopts impulsive reverse gravity 
acceleration with a low frequency of 0. 1 Hz for the activation of fluid 
reorientation with liquid filled levels of 30, 50, 65, 70 and 80«. 

II. Numerical Simulation of Liquid Hydrogen 
Reorientation with Geyser Initiation at Low Frequency 
Impulsive Reverse Gravity Acceleration of 0.1 Hz 


The present study examines time-dependent fluid behaviors, in 
particular the dynamics of liquid hydrogen and hydrogen vapor 
reorientation induced by reverse gravity acceleration which is great 
enough to introduce geyser initiation. As in Paper I 9 , time- 
dependent axial symmetry mathematical formulation are adopted. 
Detailed description of mathematical formulation, initial and 
boundary conditions suitable for the analysis of cryogenic fluid 
management under microgravity environment are given in our earlier 


80 


studies. 4,6,1 1-13 The initial profiles of liquid-vapor interface 
are determined from computations based on algorithms developed for 
the steady state formultion of microgravity fluid management. 3-7 

Detailed description of computational algorithm applicable to 
microgravity fluid management are illustrated in Paper I 9 and our 
earlier studies. 4 • 6 • 1 1 " 1 3 As we have indicated in Paper I 9 , for the 
purpose of facilitating easy comparison between computational 
resuslts and experimental measurement, a model of 0.01 size prototype 
is adopted in the computer simulation. Model size is height L - 
4.23672 cm and diameter D = 4.2672 cm. If the spacecraft had been 
coasting for a long time, aligned with its direction of motion, the 
most signigicant force, drag, would be axial and with acceleration of 
10 -4 g 0 along upward direction. The hydrogen vapor is, thus, 

v 

originally positioned at the bottom of the tank. The requirement to 
settle or to position liquid fuel over the outlet end of the spacecraft 
propellant tank prior to main engine restart poses a microgravity 
fluid behavior problem. 11 Retromaneuvers of spacecraft, such as 
STV, require settling or reorientation of the propellant prior to main 
engine firing. 1,11 Cryogenic liquid propellant shall be positioned 
over the tank outlet by using auxiliary thrusters (or idle-mode 
thrusters from the main engine) which provide a thrust parallel to the 
tank's major axis in the direction of flight. Similar to Paper I 9 , a 
small value of reverse gravity acceleration (downward direction) is 
provided by the propulsive thrust of small auxiliary engine to 
initiate the reorientation of liquid propellant. This small value of 
reverse gravity acceleration of propulsive thrust increases 
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gradually till reaching the critical value on which initiation of 
geyser is detected during the time period of fluid resettlement. We 
term this reverse gravity acceleration of propulsive thrust, which is 
capable to initiate geyser, as "geyser initiation gravity-level". 
Dliis geyser initiation gravity level has been investigated through 
the method of trial and error for the various liquid-filled levels as a 
base to simulate impulsive reverse gravity acceleration with 
frequencies of 0.1, 1.0 and 10 Hz. As we have indicated in Paper I 9 , 
cryogenic liquid hydrogen at temperature of 20K is considered. 
Hydrogen density of 0.071 g/cm 3 ; surface tension coefficient at the 
interface between liquid hydrogen and hydrogen vapor of 1.9 dyne/ cm, 
hydrogen viscosity coefficient of 1.873 x 10 3 cm 2 /s; and contact 
angle of 0.5° are used in the computer simulation. 

In this paper, among three categories of impulsive reverse 
gravity acceleration with frequencies of 0.1, 1.0 and 10 Hz, 
reorientation of cryogenic liquid hydrogen activated by geyser 
initiation impulsive reverse gravity acceleration with a low 
frequency of 0 . 1 Hz produced by propulsive thrust will be investigated 
for various liquid filled levels of propellant tank. Paper I shows 
that these geyser initiation gravity levels are 5.5 x 10 2 , 6.52 x 
10“ 2 , 6.6 x 10“ 2 , 6.7 x 10“ 2 and 8.2 x 10‘ 2 g 0 for liquid filled levels 
of 30, 50, 65, 70, and 80%, respectively. 

Table 1 shows some basic geometrys and characteristics of 
cryogenic liquid hydrogen resettlement activated by reverse gravity 
acceleration at geyser initiation gravity level. Average liquid 
height h, and maximum liquid height h m are shown in Figure 1. Average 
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free fall velocity V f , average free fall time t f , and free fall 
velocity from maximum liquid height V f]B are computed from the 
following equations: ' 


v f = 

(2g i0 h) 1/2 
2h \ 1/2 

(1) 

rt | 
II 

'9 i o * 

(2) 

II 

m 

> 

(29i „h „) 1/2 

(3) 


where g io denotes geyser initiation reverse gravity acceleration. 

To show examples of the selected sequences of time evolution of 
fluid reorientation for cryogenic hydrogen, Figures 2, 3, 4, 5, and 6 
show time evolution of fluid reorientation activated by geyser 
initiation impulse reverse gravity acceleration with a frequency of 
0.1 Hz for liquid filled levels of 30, 50, 65, 70 and 80 %, 
respectively. Each figure contains four sub-figures. Subfigure 
(A) is initial profile of liquid-vapor interface at the moment of the 
starting of fluid reorientation at time t = 0; subfigure (B) , the flow 
profile during the course of fluid reorientation before the 
initiation of geysering motion; subfigure (C) , the flow profile with 
geysering motion; and subfigure (D) , the flow profile after the ending 
of geysering motion. 

Examples of selected sequences of time evolution of fluid 
reorientation illustrate following flow behaviors: (1) The liquid 
starts to flow in an annular sheet along the solid wall of tank and 
gradually pushes the vapor toward the central portion of the lower 
dome of tank as the net acceleration, reversing the direction of 
gravity field, which is applied toward the downward direction of the 
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tank's major axis, by using small auxiliary thrusters; (2) As the 
downward fluid annular sheet along the tank wall reaches the central 
bottom dome side of the tank, a geysering flow is observed; and (3) The 
vapor is thus pushed upward centrally into the liquid and the 

geysering disappears . 

Based on the computer simulation of flow field values of maximum 
flow velocity V m , time for observing maximum flow velocity t m , and 
time for reorienting liquid flowing down and reaching the bottom of 
propellant tank t R , are obtained and illustrated in Table 2 for 
reverse gravity acceleration with impulsive frequency of 0.1 Hz. 
Scale length of maximum flow velocity L m , and scale flow acceleration 
associated with maximum velocity a m , can be computed from the 
following parameters: 




Following dimensionless parameters are introduced: V m /V f , t R /t f , 
Vt„ a m /a g , L m /h and V m /V fm where a g stands geyser initiation 
acceleration (cm/s 2 ) for corresponding geyser initiation gravity 
level g i0 . Impulsive reverse gravity acceleration, g { with 
frequency f Hz is defined as follows: 


g L = <?i 


1 + — sin 27tft 
2 


( 6 ) 


Figures 7 to 9 show the variations of dimensionless parameters in 
terms of liquid filled levels for impulsive acceleration with 
frequency of 0.1 Hz. Denominators of these six dimensionless 
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parameters are either predetermined from the geometry of liquid fill 
levels or can be deduced from the corresponding calculations 
associated with the geyser initiation gravity levels. 
Characteristics of these near constant values dimensionless 
parameters can provide a good understanding of the physics of 
microgravity fluid behaviors, in particular the active category of 
liquid acquisition or positioning, and also the design criteria of on- 
orbit spacecraft propulsion system at the critical value of reverse 
gravity acceleration of propulsive thrust which is capable to 

initiate geyser. 

An efficient propellant settling technique should minimize 
propellant usage and weight penalties. This can be accomplished by 
providing optimal acceleration to the spacecraft such that the 
propellant is reoriented over the tank outlet without any vapor 
entrainment, any excessive geysering, or any other undesirable fluid 
motion. in particular, it is important to study how well the 
impulsive acceleration can provide higher efficient propellant 
settling technique than the constant acceleration thrust technique. 
Also, what is the most optimal choice of impulsive frequency which can 
achieve the best fluid acquisition management being the goal of our 

research. 

Figure 7(A) shows the ratio of maximum flow velocity to average 
free fall flow velocity V ra /V f and its associated parameters of V m and 
V f in terms of liquid filled levels for impulsive acceleration with 
frequency of 0.1 Hz. It shows that the ratio of V m /V f varies in the 
range of 5 . 0 to 5 . 1 in the entire liquid filled levels while V m andV f 
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vary from 64.0 to 85.3 cm/s (decreasing with increasing liquid filled 
levels) and from 12.5 to 17.0 cm/s (also decreasing with increasing 
liquid filled levels) respectively . As V f can be predetermined from 
geyser initiation gravity level and average liquid height, shown in 
Equation (1) , one can make an approximate prediction of maximum flow 
velocity during the liquid reorientation for the various liquid 
filled levels. In comparison between impulsive acceleration with 
frequency of 0.1 Hz and constant thrust acceleration, it shows that 
impulsive acceleration can produce higher maximum flow velocity 

(V B /V f = 5.0 to 5.1) than that of constant thrust acceleration (V m /V f - 

4.3 to 4.9). 

Figure 7(B) shows the ratio of liquid reaching bottom time to 
average free fall time t R /t f and its associated parameters of t R and t f 
in terms of liquid filled levels for impulsive acceleration with 
frequency of 0.1 Hz. It shows that the ratio of t R /t f varies in the 
range of 1.1 to 1.2 in the entire liquid filled levels while t R and t f 
vary from 0.19 to 0.3 6 s (decreasing with increasing liquid filled 
levels) and from 0.15 to 0.32 s (also decreasing with increasing 
liquid filled levels) , respectively. As t f can be predetermined from 
geyser initiation gravity level and average liquid height, shown in 
Equation (2), one can predict the time reorienting liquid fluid 
flowing down from the original position and reaching the bottom of 
propellant tank for the various liquid filled levels at the reverse 
gravity acceleration capable for the initiation of geyser. In 
comparison between impulsive acceleration with frequency of 0.1 Hz 
and constant thrust acceleration, it shows that impulsive 


acceleration activate the flow which takes shorter time for flow to 

reach the tank bottom (t R /t f - 1.1 to 1.2) than that of constant thrust 

acceleration (t R /t f = 1.21 to 1.30). 

Figure 8 (A) shows the ratio of time for observing maximum flow 
velocity to average free fall time t./t f and its associated parameters 
of t. and t f in terms of liquid filled levels for impulsive 
acceleration with frequency of 0.1 Hz. It shows that the ratio of 
t /L varies in the range of 1.2 to 1.3 in the entire liquid filled 
levels while t„ and t f vary from 0.18 to 0.39 s (decreasing with 
increasing liquid filled levels) and from 0.15 to 0.32 s (also 
decreasing with increasing liquid filled levels), resectively. As 
we indicated in Figure 7(B) , t f can be predetermined, one can predict 
the time for observing maximum flow velocity for various liquid filled 
levels at the reverse gravity acceleration capable for the initiation 
of geyser. In comparison between impulsive acceleration with 
frequency of 0.1 Hz and constant thrust acceleration, it shows that 
impulsive acceleration can activate higher maximum flow velocity 
within a shorter period of time ( 0.18 to 0.39 s) than that of constant 

thrust acceleration (0.20 to 0.425). 

Figure 8(B) shows the ratio of scale flow acceleration 

associated with maximum velocity to geyser initiation acceleration 
for corresponding geyser initiation gravity level a m /a g and its 
associated parameters of a m and a g in terms of liquid filled levels for 
impulsive acceleration with frequency of 0.1 Hz. It shows that the 
ratio of a m /a g varies in the range of 4 . 1 to 4 . 4 in the entire liquid 
filled levels while a ra and a g vary from 218 to 355 cm/s* (increasing 
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with increasing liquid filled levels) and from 53.9 to 80.4 cm/s 2 
(also increasing with increasing liquid filled levels) , 
respectively. As a g -'can be predetermined from geyser initiation 
gravity level, one can make an approximate prediction of scale flow 
acceleration associated with maximum velocity , which is defined in 
Equation (5) , at the reverse gravity acceleration capable for the 
initiation of geyser. In comparison between impulsive acceleration 
with frequency of 0.1 Hz and constant thrust acceleration, it shows 
that impulsive acceleration can produce higher scale flow 
acceleration associated with maximum flow velocity (a m /a g = 4.1 to 
4.4) than that of constant thrust acceleration (a m /a g =3.3 to 3.8) . 

Figure 9 (A) shows the ratio of scale length of maximum flow 
velocity to average liquid height L m /h and its associated parameters 
of L ra and h in terms of liquid filled levels for impulsive acceleration 
with frequency of 0.1 Hz. It shows that the ratio of L ra /h varies in 
the range of 12.3 to 12.6 in the entire liquid filled levels while L m 
and h vary from 11.5 to 33.2 cm (decreasing with increasing liquid 
filled levels) and from 0.93 to 2.67 cm (also decreasing with 
increasing liquid filled levels) , respectively. As h can be 
predetermined from the geometry of liquid filled levels, one can make 
an approximate prediction of scale length of maximum flow velocity , 
which is defined in Equation (4) , at the reverse gravity acceleration 
capable for the initiation of geyser. In comparison between 
impulsive acceleration with frequency of 0.1 Hz and constant thrust 
acceleration, it shows that impulsive acceleration can produce longer 
scale length of maximum flow velocity ( L m = 11.5 to 33.2 cm) than that 


88 



of constant thrust acceleration (L m - 12.4 to 30.9). 

Figure 9 (B) shows the ratio of maximum flow velocity to free fall 
velocity from maximum liquid height V n /V f]> and its associated 
parameters of V m and V fm in terms of liquid filled levels for impulsive 
acceleration with frequency of 0.1 Hz. It shows that the ratio of 
V /V, varies in the range of 4.0 to 4.4 in the entire liquid filled 
levels while V m and V fm vary from 64.0 to 85.3 cm/s (decreasing with 
increasing liquid filled levels) and from 16.4 to 19.2 cm/s (also 
decreasing with increasing liquid filled levels) , respectively. As 
V fm can be predetermined from geyser initiation gravity level and 
maximum liquid height, shown in Equation (3), one can make an 
approximate prediction of maximum flow velocity at the reverse 
gravity acceleration capable for the initiation of geyser. In 
comparison between impulsive acceleration with frequency of 0.1 Hz 
and constant thrust acceleration, it shows that impulsive 
acceleration can produce higher maximum velocity (V m /V fin = 4 . 0 to 4 . 4) 
than that of constant thrust acceleration (V m /V fm = 3.8 to 4.0). 

As we have illustrated in Paper I 9 , six dimensionless parameters 
presented in this study show that the parameters hold near constant 
values through the entire ranges of liquid filled levels during the 
course of reorientation of liquid hydrogen activated by the reverse 
gravity acceleration which is great enough to initiate geysering 
flow. As the purpose of present study is to investigate an efficient 
propellant settling technique which is able to minimize propellant 
usage and weight penalties, comparison of flow parameters between 
impulsive acceleration with frequency of 0.1 Hz and that of constant 
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thrust acceleration have been made toward this objective. It shows 
that the operation of auxiliary engine with impulsive acceleration of 
0 . 1 Hz frequency is better than that of constant thrust acceleration 
in terms of introducing a higher maximum flow velocity, a shorter time 
period to reach maximum flow velocity, a shorter time period for flow 
to reach the tank bottom (tank outlet) for resettling, a longer scale 
length of maximum flow velocity, and a higher flow acceleration 
associated with maximum flow velocity. 

IV. Discussion and Conclusions 

The requirement to settle or to position liquid fluid over the 
outlet end of the spacecraft propellant tank prior to main engine 
restart poses a microgravity fluid behavior problem. Retromaneuvers 
of spacecraft require settling or reorientation of the propellant 
prior to main engine firing. Cryogenic liquid propellant is 
positioned over the tank outlet by using small auxiliary thrusters (or 
idle-mode thrusters from the main engine) which provide a thrust 
parallel to the tank's major axis in the direction of flight. 

The results of the study of fluid reorientation have to be 
evaluated in terms of how well they can be managed efficiently. An 
efficient propellant settling technique should minimize propellant 
usage and weight penalties through the operation of small thrusters 
(or idle-mode thrusters from the main engine) . This can be 
accomplished by providing optimal acceleration to the spacecraft such 
that the propellant is reoriented over the tank outlet without any 
vapor entrainment, any excessive geysering, or any other undesirable 
fluid motion. 
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Production of geyser during the propellant reorientation is not 
a desirable motion for the space fluid management under microgravity 
environment. It is because geyser is always accompanied by the vapor 
entrainment and globule formation. Geyser is observed at reverse 
gravity acceleration greater than certain critical values of 
acceleration during the course of liquid reorientation. In this 
paper, numerical simulation of positive liquid acquisition is 
attempted by introducing reverse gravity acceleration, resulting 
from the propulsive thrust with high impulsive frquency of 0.1 Hz 
auxiliary engine, which exceeds critical value for the initiation of 

geyser. 

Evaluation of performance is based on how efficient the 
impulsive&a+60Hreverse gravity with 0.1 Hz frequency in comparison wit 
constant thrust acceleration can activate following flow parameters 
at the same background thrust accelerations: (A) a higher maximum 

flow velocity, (B) a shorter time period for flow to reach maximum 
velocity, (C) a shorter time period for flow to reach tank bottom (tank 
outlet) for fluid resettling, (D) a larger length scale of maximum 
flow velocity, and (E) a higher flow acceleration associated with 
maximum flow velocity. Comparison between the results of present 
study for impulsive thrust with 0.1 Hz frequency and that of constant 
thrust, shown in Paper I 9 , it shows that impulsive thrust, is superior 
than the constant reverse thrust in terms of efficient operation of 
fluid reorientation. 

Based on the computer simulation of flow fields during the course 
of fluid reorientation, six dimensionless parameters are presented 
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both in Paper I 9 and this study. It is shown that these parameters 
hold near constant values through the entire ranges of liquid filled 
levels during the course of fluid reorientation activated by the 
reverse gravity acceleration in both constant and impulsive thrusts 
great enough to initiate geyser. As the denominators of these 
dimensionless parameters are either predetermined from the geometry 
of liquid filled levels, as shown in Table 1, or can be deduced from the 
corresponding calculations associated with the geyser initiation 
gravity levels, one can predict the values of these flow parameters. 
Present study can greatly enhance our understanding in the behaviors 
of cryogenic fluid resettlement under reduced gravity environment. 
This is particularly important for liquid acquisition technique to be 
used in on-orbit spacecraft design. 

Any fluid capable of motion relative to the spacecraft will be 
subject to an acceleration relative to the mass center of the 
spacecraft that arises from the gravity gradient of the Earth 1 ' 5 . 
In addition to the Earth's gravitational force, the interaction 
between the particle mass of fluids and the spacecraft mass due to 
gravity gradient acceleration 1 9 have also been taken into 
consideration in this microgravity fluid management study. 

To conclude, we have demonstrated that, the computer algorithm 
presented, can be used to simulate fluid behavior in a microgravity 
environment, in particular the development of technology necessary 
for acquisition or positioning of liquid and vapor within a tank to 
enable liquid outflow or vapor venting through active liquid 
acquisition by the creation of a positive acceleration environment 
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resulting from propulsive thrust. Better understanding of the full 
pictures of flow fields in both constant and impulsive thrusts, during 
the course of fluid 'reorientation can provide the proper design 
techniques for handling and managing the cryogenic liquid propellants 
to be used in on-orbit spacecraft propulsion. It is important to 
emphasize that impulsive reverse gravity thrust is superior to 
constant reverse gravity thrust for the activation of liquid 
necessary for the resettlement of liquid in a reduced gravity 

environement . 
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TABLE 1 SOME BASIC GEOMETRYS AND CHARACTERISTICS 
OF CRYOGENIC LIQUID HYDROGEN REORIENTATION 



(cm/s) 

































Table 2 CHARACTERISTICS OF CRYOGENIC HYDROGEN REORIENTATION 
(Frequency of Impulsive Acceleration “0.1 Hz) 






































Figure Captions 

Figure 1. (A) Distribution of grid points in the radial-axial plane 

of cylindrical coordinate for propellant tank, and (B) 
Model size propellant tank adopted for numerical 
simulation with geometrical description. 

Figure 2 . Selected sequences of time evolution of fluid 

reorientation with liquid filled level of 30% for 
impulsive acceleration with frequency of 0.1 Hz, (A) 
initial profile, (B) flow profile before the initiation of 
geyser, (C) flow profile with geyser, and (D) flow profile 
after the ending of geyser. 

Figure 3 . Selected sequences of time evolution of fluid 

reorientation with liquid filled level of 50% for 

impulsive acceleration with frequency of 0.1 Hz, (A) 

initial profile, (B) flow profile before the initiation of 
geyser, (C) flow prof ile with geyser, and (D) flow profile 
after the ending of geyser. 

Figure 4. Selected sequences of time evolution of fluid 

reorientation with liquid filled level of 65% for 

impulsive acceleration with frequency of 0.1 Hz, (A) 
initial profile, (B) flow profile before the initiation of 
geyser, (C) flow prof ile with geyser , and (D) flow profile 
after the ending of geyser. 

Figure 5. Selected sequences of time evolution of fluid 

reorientation with liquid filled level of 70% for 

impulsive acceleration with frequency of 0.1 Hz, (A) 
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Figure 6 


Figure 7 . 


Figure 8 . 


Figure 9 . 


initial profile, (B) flow profile before the initiation of 
geyser, (C) flow prof ile with geyser, and (D) flow profile 
after the 'ending of geyser. 

. Selected sequences of time evolution of fluid 
reorientation with liquid filled levels of 80% for 
impulsive acceleration with frequency of 0.1 Hz, (A) 
initial profile, (B) flow profile before the initiation of 
geyser, (C) flow prof ile with geyser, and (D) flow profile 
after the ending of geyser. 

(A) Ratio of V m /V f and its associated parameters in terms 
of liquid filled levels, (B) Ratio of t R /t f and its 
associated parameters in terms of liquid filled levels. 
(A) Ratio of t m /t f and its associated parameters in terms 
of liquid filled levels, (B) Ratio of a m /a^ and its 
associated parameters in terms of liquid filled levels. 
(A) Ratio of L m /h and its associated parameters in terms of 
liquid filled levels, (B) Ratio of V m /V fm and its 
associated parameters in terms of liquid filled levels. 
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ABSTRACT 


The requirement to settle or to position liquid fluid over the 
outlet end of spacecraft propellant tank prior to main engine restart 
poses a microgravity fluid behavior problem. Resettlement or 
reorintation of liquid propellant can be accomplished by providing 
optimal acceleration to the spacecraft such that the propellant is 
reorientaed over the tank oulet without any vapor entrainment, any 
excessive geyser ing, or any other undesirable fluid motion for the 
space fluid management under microgravity environment. The purpose 
of present study is to investigate most efficient technique for 
propellant resettling through the minimization of propellant usage 
and weight penalties. Comparison between the constant reverse 
gravtiy acceleration and impulsive reverse gravity acceleration to be 
used for the activation of propellant resettlement, it shows that 
impulsive reverse gravity thrust is superior to constant reverse 
gravity thrust for liquid reorientation in a reduced gravity 
environment. Comparison among impulsive reverse gravity thrust with 
0.1, 1.0 and 10 Hz frequencies for liquid filled level in the range 
between 30 to 80 %, it shows that the selection of 1.0 Hz frequency 
impulsive thrust over the other frequency ranges of impulsive thrust 
is most proper based on the present study. 
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Nomenclature 


a g - geyser initiation acceleration (cm/s ) 

a = scale flow acceleration associated with maximum velocity 

in 

(cm/s 2 ), defined by Equation (5) 

D = diameter of propellant tank (cm) 
f = frequency of impulsive thrust (H z ) 

g. = impulsive reverse gravity acceleration, defined by Equation (6) 
g io = geyser initiation gravity-level (g 0 ) 

2 

g 0 = normal Earth gravitational acceleration = 9.81 m/s 
h = average liquid height (cm) 

h m = maximum liquid height (cm) 

L = height of propellant tank (cm) 

L m = scale length of maximum liquid height (cm/s) , defined by 
Equation (4) 

STV = Space Transfer Vehicle 
t f = average free fall time (s) 

t m = time for observing maximum flow velocity (s) 
t R = liquid reaching tank bottom time (s) 

V fm = free fall velocity from maximum liquid height (cm/s) , defined by 
Equation (3) 

V f = average free fall velocity (cm/s), defined by Equation (2) 
V ra = maximum flow velocity (cm/s) 
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I . Introduction 


In spacecraft design, the requirements fora settled propellant 
are different for tank pressurization, engine restart, venting, or 
propellant transfer. Prepressurization requires that heat and mass 
transfer effects be minimized; otherwise, a process of chill down of 
tank, venting of noncondensing gases, etc. , may have to carry out for 
the cryogenic system. For engine restart, it is necessary to have the 
liquid settle with no bubbles near the tank outlet so that the initial 
flow of propellant will not carry vapor to the pump or engine. The 
slosh wave amplitude should be relatively low to keep the center of 
mass shifts within an acceptable range and wave motion low enough to 
avoid pressure collapse caused by interface agitation . For venting, 
it is probably necessary that virtually all bubbles be displaced from 
the bulk liquid so that a two-phase mixture is not vented. Propellant 
transfer requires that the liquid be completely settled with 
virtually no bubbles. Outflow of a liquid near the tank outlet can 
result in the prematrue ingestion of gas while a signigicant amount of 
liquid is still in the tank under microgravity environment. This 
phenomenon is termed "suction dip". Slosh wave motion must be 
minimal because the combination of "suction dip" and sloshing could 
cause gas pull-through to occur more readily in microgravity than if 

the surface were essentially quiescent. 

During the prepressurization of a cryogenic propellant in 
microgravity, significant heat and mass transfer will occur if the 
liquid interface is disturbed. Interface disturbances may result 
from (a) impingement of the gas on the liquid surface at a mass flow 


102 



rate sufficient to cause Kelvin-Helmholtz instability, (b) globule 
formation from breaking waves caused by wave motion over baffles or 
internal hardware, (cj globule and surface froth formation resulting 
from movement of bubbles through the liquid to the surface, and (d) 
surface froth formation because of gas impingement. 

Bubble and globule formation as a result of liquid impact with 
the aft end of the tank could lead to propellant loss for the 
spacecraft during venting. Globules could be entrained in the vented 
ullage gas or bubbles rising through the liquid and expanding because 
of the decreasing tank pressure could cause a spray of globules to be 
vented. Liquid level rise, vent liquid loss, fluid freezing, and 
vehicle dynamics are all affected by the microgravity levels. 

A key objective of the cryogenic fluid management of spacecraft 
propulsion system, such as a Space Transfer Vehicle 1 (STV) , is to 
develop the technology necessary for acquisition or positioning of 
liquid and vapor within a tank in reduced gravity to enable liquid 
outflow or vapor venting. Liquid acquisition techniques can be 
divided into two general categories: (1) Active liquid acquisition 
by the creation of a positive acceleration environment resulting from 
the propulsive thrust of small auxiliary engines, and (2) Passive 
liquid acquisition utilizing the liquid capillary forces provided by 
using solid baffles of liquid traps made of fine mesh screen material . 
In this series of study, active liquid acquisition is aimed for 
numerically simulating the resettlement of cryogenic liquid 
hydrogen. Liquid hydrogen, which, in general, poses more severe 
technical challenges than liquid oxygen, is used as the test bed 
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working fluid in this study. 

Recently Lesie 2 was able to measure and to numerically compute 
the bubble shapes, at various ratios of centrifugal force to surface 
tension force in 2 , 4 and 6.3 cm deep cylinders in the microgravity 
environment. The results showed excellent agreement between model 
computation and measurements. Hung and Leslie 3 extended Leslie's 
work 2 to rotating free surfaces influenced by gravity with higher 
rotating speeds when the bubble intersects with both the top and 
bottom walls of the cylinder. Hung et al. 4 ' 5 further extended the 
work to include rotating speeds which resulted with bubbles 
intersecting and/or without intersecting the top, bottom and side 
walls of the cylinder. 

An analysis of time-dependent dynamical behavior of surface 
tension on partially-filled rotating fluids in both low gravity and 
micr ogr av i ty environments was carried out by numerically solving the 
Navier-Stokes equations subjected to the initial and the boundary 
conditions 4 ' 6 . At the interface between the liquid and the gaseous 
fluids, both the kinematic surface boundary condition, and the 
interface stress conditions for components tangential and normal to 
the interface, were applied. The initial condition for the bubble 
profiles was adopted from steady-state formulations developed by Hung 
and Leslie 3 , and Hung et al . 6 for rotating cylinder tank; and by Hung 
et al. 5 ' 7 for the dewar-shaped container to be used in the Gravity 
Probe-B Spacecraft 8 . Some of the steady-state formulations of 
bubble shapes, in particular for bubbles intersecting at the top wall 
of the cylinder, were compared with the experiment carried out by 
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Leslie 2 in a free-fallinng aircraft (KC-135) . Comparisons of time- 
dependent results between numerical computations and experiments 
were unavailable. This was because the calibration of the recordings 
of time-dependent gravity variations in a KC-135 aircraft during the 
short time periods of microgravity environment is very difficult. 
There was also an unavailability of accelerometer data for measuring 
the actual levels of microgravity during the experiment. 

An efficient propellant settling technique should minimize 
propellant usage and weight penalties. This can be accomplished by 
providing optimal acceleration to the spacecraft such that the 
propellant is reoriented over the tank outlet without any vapor 
entrainment, any excessive geysering, or any other undesirable fluid 
motion . 

Production of geyser during the propellant reorientation is not 
a desirable motion for the space fluid management. It is because 
geyser is always accompanied by the vapor entrainment and globule 
formation. Geyser is observed at reverse gravity thrust greater than 
certain critical values of acceleration during the course of liquid 
reorientation. In other words, geyser will not be observed at very 
low reverse gravity level, and it will be detected when the reverse 
gravity level is greater than the certain critical value. In this 
series of study, numerical simulation of positive liquid acquisition 
is attemted by introducing reverse gravity acceleration, resulting 
from the propulsive thrust of small auxiliary engines which exceeds 
the critical value for geyser initiation. The reverse gravity 
acceleration is starting with a small value and increases gradually 
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tin the initiation of geyser is detected in the computer simulation 
f o£- the liguid reorientation of propellant hank with various liguid - 
filled levels. 

In this series of studies , time-dependent computations have been 
carried out to investigate the dynamical behavior of fluid 
reorientation or resettling of propellant prior to main engine firing 
for spacecraft restart at the net reverse gravity acceleration which 
is great enough to initiate geyser during the liquid reorientation. 
First paper of present study (Paper I ) 9 investigates the 
characteristics of fluid resettlement due to the reversal of lowest 
constant reverse gravity acceleration which is great enough to 
initiate geyser without any impulsive acceleration. The frequency 
of impulsive acceleration is generally termed "gravity jitters". 
Gravity jitters are produced by spacecraft attitude motion, machinery 
(turbine, pump, engine) vibrations, thruster firing, thruster 
shutdown, impulsive engine acceleration, etc . 10 Positioning of 
liquid propellant over the tank outlet can be carried out by using 
small auxiliary thrusters which provide a thrust parallel to the 
tank' s major axis in the direction of flight. 

Computer simulation of flow field based on Paper I during the 
course of fluid reorientation induced by constant reverse gravity 
acceleration show that six dimensionless parameters are abtained in 
the study. These parameters hold near constant values through the 
entire ranges of liquid filled levels during the course of fluid 
reorientation activated by the reverse gravity acceleration great 
enough to initiate geyser. As the denominators of these 
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dimensionless parameters are either predetermined from the geometry 
of liquid filled levels or can be deduced from the corresponding 
calculations associated with the geyser initiation gravity levels. 
One can predict the values of these flow parameters. These 
predictable parameters include maximum flow velocity V m , time for 
observing maximum flow velocity t m , time for reorienting liquid 
flowing down and reaching the bottom of propellant tank t R/ scale 
length of maximum flow velocity L B , and scale flow acceleration 

associated with maximum velocity a B . 

Instead of applying constant reverse gravity acceleration as we 
described in Paper I 9 , this paper adopts impulsive reverse gravity 
acceleration with a medium frequency of 1.0 Hz for the activation of 
fluid reorientation with liquid filled levels of 30, 50, 65, 70 and 

80%. 


IX. Numerical Simulation of Liquid Hydrogen 
Reorientation with Geyser Initiation at Medium 
Frequency Impulsive Reverse Gravity Acceleration of 1.0 Hz 

The present study examines time-dependent fluid behaviors, in 
particular the dynamics of liquid hydrogen and hydrogen vapor 
reorientation induced by reverse gravity acceleration which is great 
enough to introduce geyser initiation. As in Paper I 9 , time- 
dependent axial symmetry mathematical formulation are adopted. 
Detailed description of mathematical formulation, initial and 
boundary conditions suitable for the analysis of cryogenic fluid 
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management under microgravity environment are given in our earlier 
studies. 4 ' 6 ' 1 1_ 13 The initial profiles of liquid-vapor interface 
are determined from computations based on algorithms developed for 
the steady state formulation of microgravity fluid management. 3 7 

Detailed description of computational algorithm applicable to 
microgravity fluid management are illustrated in Paper X and our 
ier studies. As we have indicated in Paper X , for the 

purpose of facilitating easy comparison between computational 
resuslts and experimental measurement, a model of 0.01 size prototype 
is adopted in the computer simulation. Model size is height L = 
4.23672 cm and diameter D = 4.2672 cm. If the spacecraft had been 
coasting for a long time, aligned with its direction of motion, the 
most signigicant force, drag, would be axial and with acceleration of 
10“ 4 g 0 along upward direction. The hydrogen vapor is, thus, 
originally positioned at the bottom of the tank. The requirement to 
settle or to position liquid fuel over the outlet end of the spacecraft 
propellant tank prior to main engine restart poses a microgravity 
fluid behavior problem. 11 Retromaneuvers of spacecraft, such as 
STV, require settling or reorientation of the propellant prior to main 
engine firing. 1/11 Cryogenic liquid propellant shall be positioned 
over the tank outlet by using auxiliary thrusters (or idle-mode 
thrusters from the main engine) which provide a thrust parallel to the 
tank's major axis in the direction of flight. Similar to Paper I , a 
small value of reverse gravity acceleration (downward direction) is 
provided by the propulsive thrust of small auxiliary engine to 
initiate the reorientation of liquid propellant. This small value of 
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reverse gravity acceleration of propulsive thrust increases 
gradually till reaching the critical value on which initiation of 
geyser is detected during the time period of fluid resettlement- We 
term this reverse gravity acceleration of propulsive thrust, which is 
capable to initiate geyser, as "geyser initiation gravity-level". 
This geyser initiation gravity level has been investigated through 
the method of trial and error for the various liquid-filled levels as a 
base to simulate impulsive reverse gravity acceleration with 
frequencies of 0.1, 1.0 and 10 Hz. As we have indicated in Paper I 9 , 
cryogenic liquid hydrogen at temperature of 20K is considered. 
Hydrogen density of 0.071 g/cm 3 ; surface tension coefficient at the 
interface between liquid hydrogen and hydrogen vapor of 1.9 dyne/ cm; 
hydrogen viscosity coefficient of 1.873 x 10" 3 cm 2 /s; and contact 
angle of 0.5° are used in the computer simulation. 

In this paper, among three categories of impulsive reverse 
gravity acceleration with frequencies of 0.1, 1.0 and 10 Hz, 
reorientation of cryogenic liquid hydrogen activated by geyser 
initiation impulsive reverse gravity acceleration with a frequency of 
1.0 Hz produced by propulsive thrust will be investigated for various 
liquid filled levels of propellant tank. Paper I shows that these 
geyser initiation gravity levels are 5.5 x 10" , 6.52 x 10 , 6.6 x 
10“ 2 , 6.7 x 10“ 2 and 8.2 x 10" 2 g 0 for liquid filled levels of 30, 50, 
65, 70, and 80%, respectively. 

Table 1 shows some basic geometrys and characteristics of 
cryogenic liquid hydrogen resettlement activated by reverse gravity 
acceleration at geyser initiation gravity level. Average liquid 
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height h, and maximum liquid height h m are shown in Figure 1. Average 
free fall velocity V f , average free fall time t f , and free fall 
velocity from maximum liquid height V fm are computed from the 
following equations: 


(2g i 0 h ) 1 1 2 

(1) 

'2h ' 

1 / 2 

(2) 

Ui<J 



(2g i( 

h ) 1 / 2 
s AA Tn t 

(3) 


where g io denotes geyser initiation reverse gravity acceleration. 

To show examples of the selected sequences of time evolution of 
fluid reorientation for cryogenic hydrogen. Figures 2, 3, 4, 5, and 6 
show time evolution of fluid reorientation activated by geyser 
initiation impulse reverse gravity acceleration with a frequency of 
1.0 Hz for liquid filled levels of 30, 50, 65, 70 and 80 %, 

respectively. Each figure contains four sub-figures. Subfigure 
(A) is initial profile of liquid-vapor interface at the moment of the 
starting of fluid reorientation at time t = 0; subfigure (B) , the flow 
profile during the course of fluid reorientation before the 
initiation of geysering motion; subfigure (C) , the flow profile with 
geysering motion; and subfigure (D) , the flow profile after the ending 
of geysering motion. 

Examples of selected sequences of time evolution of fluid 
reorientation illustrate following flow behaviors: (1) The liquid 

starts to flow in an annular sheet along the solid wall of tank and 
gradually pushes the vapor toward the central portion of the lower 
dome of tank as the net acceleration, reversing the direction of 
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gravity field, which is applied toward the downward direction of the 
tank' s major axis, by using small auxiliary thrusters; (2) As the 
downward fluid annular sheet along the tank wall reaches the central 
bottom dome side of the tank, a geysering flow is observed; and (3) The 
vapor is thus pushed upward centrally into the liquid and the 
geysering disappears. 

Based on the computer simulation of flow field values of maximum 
flow velocity V m , time for observing maximum flow velocity t m , and 
time for reorienting liquid flowing down and reaching the bottom of 
propellant tank t R , are obtained and illustrated in Table 2 for 
reverse gravity acceleration with impulsive frequency of 1.0 Hz. 
Scale length of maximum flow velocity L m , and scale flow acceleration 
associated with maximum velocity a m , can be computed from the 
following parameters: 




V 


m 


a m = 


(4) 

(5) 


Following dimensionless parameters are introduced: 

t ra /t f , a m /a g , L m /h and V m /V fm where a g stands geyser initiation 
acceleration (cm/s 2 ) for corresponding geyser initiation gravity 
level 9 i o • Impulsive reverse gravity acceleration, 9j with 
frequency f Hz is defined as follows: 


Si - Si 


i o 


1 ' 
1 + — sin 2irft 
2 


(6) 


Figures 7 to 9 show the variations of dimensionless parameters in 
terms of liquid filled levels for impulsive acceleration with 
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Denominators of these six dimensionless 


frequency of 1.0 Hz. 
parameters are either predetermined from the geometry of liquid fill 
levels or can be deduced from the corresponding calculations 
associated with the geyser initiation gravity levels. 
Characteristics of these near constant values dimensionless 
parameters can provide a good understanding of the physics of 
microgravity fluid behaviors, in particular the active category of 
liquid acquisition or positioning, and also the design criteria of on- 
orbit spacecraft propulsion system at the critical value of reverse 
gravity acceleration of propulsive thrust which is capable to 
initiate geyser. 

An efficient propellant settling technique should minimize 
propellant usage and weight penalties. This can be accomplished by 
providing optimal acceleration to the spacecraft such that the 
propellant is reoriented over the tank outlet without any vapor 
entrainment, any excessive geysering, or any other undesirable fluid 
motion. In particular, it is important to study how well the 
impulsive acceleration can provide higher efficient propellant 
settling technique than the constant acceleration thrust technique. 
Also, what is the most optimal choice of impulsive frequency which can 
achieve the best fluid acquisition management being the goal of our 
research. 

Figure 7 (A) shows the ratio of maximum flow velocity to average 
free fall flow velocity V n /V f and its associated parameters of V m and 
V f in terms of liquid filled levels for impulsive acceleration with 
frequency of 1.0 Hz. It shows that the ratio of V ra /V f varies in the 
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range of 6.6 to 6.7 in the entire liquid fillsd levels while V m and Vp 
vary from 82 . 5 to 112 . 2 cm/s (decreasing with increasing liquid filled 
levels) and from 12.5' to 17.0 cm/s (also decreasing with increasing 
liquid filled levels) , respectively. As V f can be predetermined from 
geyser initiation gravity level and average liquid height, shown in 
Equation (1) , one can make an approximate prediction of maximum flow 
velocity during the liquid reorientation for the various liquid 
filled levels. In comparison between impulsive acceleration with 
frequency of 1 . 0 Hz and with frequencies of 0 . 1 and 10 Hz , it shows that 
the former case of thrust can activate higher maximum flow velocity 
(V m /V f = 6.6 to 6.7) than that of the latter two cases of thrust 
accelerations (V m /V f = 5 . 0 to 5 . 2 for f= 0.1 Hz ; and =5.1 to 5 . 3 for f — 
10 Hz) . 

Figure 7(B) shows the ratio of liquid reaching bottome time to 
average free fall time t R /tf and its associated parameters of t R and t f 
in terms of liquid filled levels for impulsive acceleration with 
frequency of 1.0 Hz. It shows that the ratio of t R /t f varies in the 
range of 1.1 to 1.2 in the entire liquid filled levels while t R and t f 
vary from 0.19 to 0.35 s (decreasing with increasing liquid filled 
levels) , respectively. As t f can be predetermined from geyser 
initiation gravity level and average liquid height, shown in Equation 
(2) , one can predict the time for reorienting liquid fluid flowing 
down from the original position and reaching the bottom of propellant 
tank for the various liquid filled levels at the reverse gravity 
acceleration capable for the initiation of geyser. In comparison 
between impulsive acceleration with frequency of 1.0 Hz and with 
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frequencies of 0.1 and 10 Hz, it shows that the former case of thrust 
can activate the flow which takes sliqhtly shorter time period for 
flow to reach the tank bottom (t R = 0.19 to 0.35 s) than that of the 
latter two cases of acceleration (t R = 0.19 to 0.36 s for f = 0 . 1 Hz ; and 

= 0.18 to 0.37 s for f =10 Hz). 

Figure 8(A) shows the ratio of time for observing maximum flow 
velocity to average free fall time t m /t f and its associated parameters 
of t n andt f in terms of liquid filled levels. It shows that the ratio 
of t m /t f varies in the range of 1.1 to 1.2 in the entire liquid filled 
levels while t m and t f vary from 0.18 to 0.35 s (decreasing with 
increasing liquid filled levels) and from 0.15 to 0.32 s (also 
decreasing with increasing liquid filled levels) , respectively. As 
we indicated in Figure 7(B) , t f can be predetermined, one can predict 
the time for observing maximum flow velocity for various liquid filled 
levels at the reverse gravity acceleration capable for the initiation 
of geyser. In comparison between impulsive acceleration with 
frequency of 1 . 0 Hz and with frequencies of 0 . 1 and 10 Hz , it shows that 
the former case of thrust can produce higher maximum flow velocity 
with a slightly shorter period of time (t m = 0. 18 to 0 . 35 s) than that 
of the latter two cases of acceleration (t B = 0.18 to 0.39 s for both f - 
0. 1 and 10 Hz) . 

Figure 8(B) shows the ratio of scale flow acceleration 
associated with maximum velocity to geyser initiation acceleration 
for corresponding geyser initiation gravity level a m /ag and its 
associated parameters of a m and a g in terms of liquid filled levels for 
impulsive aceeleration with frequency of 1.0 Hz. It shows that the 
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ratio of a„/aj varies in the range of 5.7 to 6.1 in the entire liquid 
filled levels while a m and ag vary from 320 to 458 cm/s (increasing 
with increasing liquid filled levels) and from 53.9 to 80.4 cm/s 
(also increasing with increasing liquid filled- levels), 
respectively. As a g can be predetermined from geyser initiation 
gravity level, one can make an approximate prediction of scale flow 
acceleration associated with maximum velocity, which is defined in 
Equation (5) , at the reverse gravity acceleration capable for the 
initiation of geyser. In comparison between impulsive acceleration 
with frequency ot 1 . 0 Hz and with frequencies of 0 . 1 and 10 Hz , it shows 
that the former case of thrust can activate higher scale flow 
acceleration associated with maximum flow velocity (a m /a g = 5.7 to 
6 _ i) than that of the latter two cases of acceleration (a m /a g = 4.1 to 
4.4 for f= 0.1 Hz? and = 4.2 to 4.3 for f = 10 Hz) . 

Figure (9A) shows the ratio of scale length of maximum flow 
velocity to average liquid height L m /h and its associated parameters 
of L m and h in terms of liquid levels for impulsive acceleration with 
frequency of 1.0 Hz. It shows that the ratio of L ra /h varies in the 
range of 14.7 to 15.9 in the entire liquid filled levels while L m and h 
vary from 14.8 to 39.3 cm (decreasing with increasing liquid filled 
levels) and from 0.93 to 2.67 cm (also decreasing with increasing 
liquid filled levels) , respectively. As h can be predetermined from 
the geometry of liquid filled levels, one can make an approximate 
prediction of scale length of maximum flow velocity , which is defined 
in Equation (4) , at the reverse gravity acceleration capable for the 
initiation of geyser. In comparison between impulsive acceleration 
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with frequency of 1 . 0 Hz and with frequencies of 0 . 1 and 10 Hz , it shows 
that the former case of thrust can activate larger scale length of 
maximum flow velocity *(L m /h = 14 . 7 to 15 . 9) than that of the latter two 
cases of acceleration (L m /h = 12 . 3 to 12 . 6 for f = 0 . 1 Hz and = 12 . 1 to 
13.1 for f = 10 Hz) . 

Figure 9 (B) shows the ratio of maximum flow velocity to free fall 
velocity from maximum liquid height V m /V fm and its associated 
parameters of V m and V fm in terms of liquid filled levels for impulsive 
acceleration with frequency of 1.0 Hz. It shows that the ratio of 
V /V t varies in the range of 5.2 to 5.8 in the entire liquid filled 
levels while V m and V fm vary from 82.5 to 112.2 cm/s (decreasing with 
increasing liquid filled levels) and from 16.4 to 19.2 cm/s (also 
decreasing with increasing liquid filled levels) , respectively. As 
V f m can be predetermined from geyser initiation gravity level and 
maximum liquid height, shown in Equation (3) , one can make an 
approximate prediction of maximum flow velocity at the reverse 
gravity acceleration capable for the initiation of geyser. In 
comparison between impulsive acceleration with frequency of 1.0 Hz 
and with frequencies of 0.1 and 10 Hz, it shows that the former case of 
thrust can activate higher maximum flow velocity (V m /V f m =5.2 to 5 . 8) 
than that of the latter two cases of acceleration (V m /V fm =4.0 to 4 . 4 
for f = 0.1 Hz; and = 3.9 to 4.6 for f = 10 Hz) . 

As we have mentioned in Paper I 9 , six dimensionless parameters 
presented in this study show that the parameters hold near constant 
values through the entire ranges of liquid filled levels during the 
course of reorientation of liquid hydrogen activated by the reverse 
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gravity acceleration with impulsive thrust frequency of 1.0 Hz, which 
is great enough to initiate geysering flow. As the purpose of present 
study is to investigate an efficient propellant settling technique 
which is able to minimize propellant usage and weight penalties, 
comparison of flow parameters between impulsive acceleration with 
frequency of 1.0 Hz and that of 0.1 and 10 Hzs have been made toward 
this objective. It shows that the operation of auxiliary engine with 
impulsive acceleration of 1.0 Hz frequency is better than that of 0.1 
and 10 Hz frequency impulsive thrust acceleration in terms of 
introducing a higher maximum flow velocity, a shorter time period fo 
flow to reach maximum flow velocity, a shorter time period for flow to 
reach the tank bottom (tank outlet) for fluid resettling, a larger 
scale length of maximum flow velocity, and a higher flow acceleration 
associated with maximum flow velocity. It is shown that the 
impulsive acceleration with frequencies of 0.1 and 10 Hz can provide a 
higher efficient propellant resettling technique than that of the 
constant thrust acceleration based on same arguments mentioned 
ear iier. It can be ranked in the order of efficienty that the 
impulsive acceleration with frequency of 1.0 Hz, followed by that with 
impulsive frequencies of 0. 1 and 10 Hz which are about the same order, 
and then by constant thrust acceleration based on priority of 
selecting higher efficient propellant resettling technique. 

IV. Discussion and Conclusions 
The requirement to settle or to position liquid fluid over the 
outlet end of the spacecraft propellant tank prior to main engine 
restart poses a microgravity fluid behavior problem. Retromaneuvers 
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of spacecraft require settling or reorientation of the propellant 
prior to main engine firing. Cryogenic liquid propellant is 
positioned over the tank outlet by using small auxiliary thrusters (or 
idle-mode thrusters from the main engine) which provide a thrust 
parallel to the tank's major axis in the direction of flight. 

The results of the study of fluid reorientation have to be 
evaluated in terms of how well they can be managed efficiently. An 
efficient propellant settling technique should minimize propellant 
usage and weight penalties through the operation of small thrusters 
(or idle-mode thrusters from the main engine) . This can be 
accomplished by providing optimal acceleration to the spacecraft such 
that the propellant is reoriented over the tank outlet without any 
vapor entrainment, any excessive geysering, or any other undesirable 
fluid motion. 

Production of geyser during the propellant reorientation is not 
a desirable motion for the space fluid management under microgravity 
environment. It is because geyser is always accompanied by the vapor 
entrainment and globule formation. Geyser is observed at reverse 
gravity acceleration greater than certain critical values of 
acceleration during the course of liquid reorientation. In this 
paper, numerical simulation of positive liquid acquisition is 
attempted by introducing reverse gravity acceleration, resulting 
from the propulsive thrust with impulsive frquencies of 0.1, 1.0 and 
10 Hz auxiliary engines, which exceeds critical value for the 

initiation of geyser. 

Evaluation of performance is based on how efficient the 
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impulsive reverse gravity with various frequencies can activate 
following flow parameters at the same background thrust 
accelerations: (A) a higher maximum flow velocity, (B) a shorter 
time period for flow to reach maximum velocity, (C) a shorter time 
period for flow to reach tank bottom (tank outlet) for fluid 
resettling, (D) a larger length scale of maximum flow velocity, and 
(E) a higher flow acceleration associated with maximum flow velocity . 
in summary, it is shown that impulsive reverse gravity with various 
frequencies is the most efficient choice for the reverse thrust to be 
operated at auxiliary engine for the purpose of fluid reorientation. 
Comparison between the results of present study for impulsive thrust 
with various frequencies and that of constant thrust, shown in Paper 
I 9 , it shows that impulsive thrust, regardless of frequency range, is 
always superior than the constant reverse thrust in terms of 
efficient operation of fluid reorientation. 

Based on the computer simulation of flow fields during the 
course of fluid reorientation, six dimensionless parameters are 
presented both in Paper I 9 and this study. It is shown that these 
parameters hold near constant values through the entire ranges of 
liquid filled levels during the course of fluid reorientation 
activated by the reverse gravity acceleration in both constant and 
impulsive thrusts great enough to initiate geyser. As the 
denominators of these dimensionless parameters are either 
predetermined from the geometry of liquid filled levels, as shown in 
Table 1, or can be deduced from the corresponding calculations 
associated with the geyser initiation gravity levels, one can predict 
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the values of these flow parameters. Present study can greatly 
enhance our understanding in the behaviors of cryogenic fluid 
resettlement under' reduced gravity environment. This is 
particularly important for liquid acquisition technique to be used in 
on-orbit spacecraft design. 

Any fluid capable of motion relative to the spacecraft will be 
subject to an acceleration relative to the mass center of the 
spacecraft that arises from the gravity gradient of the Earth 14 ' 15 . 
In addition to the Earth's gravitational force, the interaction 
between the particle mass of fluids and the spacecraft mass due to 
gravity gradient acceleration 1 4 have also been taken into 
consideration in this microgravity fluid management study. 

To conclude, we have demonstrated that, the computer algorithm 
presented, can be used to simulate fluid behavior in a microgravity 
environment, in particular the development of technology necessary 
for acquisition or positioning of liquid and vapor within a tank to 
enable liquid outflow or vapor venting through active liquid 
acquisition by the creation of a positive acceleration environment 
resulting from propulsive thrust. Better understanding of the full 
pictures of flow fields in both constant and impulsive thrusts, during 
the course of fluid reorientation can provide the proper design 
techniques for handling and managing the cryogenic liquid propellants 
to be used in on-orbit spacecraft propulsion. It is important to 
emphasize that impulsive reverse gravity thrust is superior to 
constant reverse gravity thrust for the activation of liquid 
necessary for the resettlement of liquid in a reduced gravity 
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environement. It is also worthwhile to mention that the selection of 
1.0 Hz frequency impulsive thrust over the other frequency ranges is 
most proper based on the present study. 
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TABLE 1 SOME BASIC GEOMETRYS AND CHARACTERISTICS 
OF CRYOGENIC LIQUID HYDROGEN REORIENTATION 



(cm/ s ) 



































Table 2 CHARACTERISTICS OF CRYOGENIC HYDROGEN REORIENTATION 
(Frequency of Impulsive Acceleration “ 1*0 Hz) 






































Figure Captions 


Figure 1. (A) Distribution of grid points in the radial-axial plane 

of cylindrical coordinate for propellant tank, and (B) 
Model size propellant tank adopted for numerical 
simulation with geometrical description. 

Figure 2. Selected sequences of time evolution of fluid 

reorientation with liquid filled level of 30% for 
impulsive acceleration with frequency of 1.0 Hz, (A) 
initial profile, (B) flow profile before the initiation of 
geyser, (C) flow prof ile with geyser, and (D) flow prof ile 
after the ending of geyser. 

Figure 3 . Selected sequences of time evolution of fluid 

reorientation with liquid filled level of 50% for 
impulsive acceleration with frequency of 1.0 Hz, (A) 
initial profile, (B) flow profile before the initiation of 
geyser, (C) flow prof ile with geyser, and (D) flow profile 
after the ending of geyser. 

Figure 4. Selected sequences of time evolution of fluid 

reorientation with liquid filled level of 65% for 
impulsive acceleration with frequency of 1.0 Hz, (A) 
initial profile, (B) flow prof ile before the initiation of 
geyser, (C) flow prof ile with geyser , and (D) flow profile 
after the ending of geyser. 

Figure 5. Selected sequences of time evolution of fluid 

reorientation with liquid filled level of 10s for 
impulsive acceleration with frequency of 1.0 Hz, (A) 
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Figure 6 


Figure 7 . 


Figure 8. 


Figure 9 . 


initial profile, (B) flow profile before the initiation of 
geyser, (C) flow prof ile with geyser, and (D) flow profile 
after the ‘ending of geyser. 

Selected sequences of time evolution of fluid 
reorientation with liquid filled levels of 80% for 
impulsive acceleration with frequency of 1.0 Hz, (A) 
initial profile, (B) flow profile before the initiation of 
geyser, (C) flow prof ile with geyser , and (D) flow prof ile 
after the ending of geyser. 

(A) Ratio of V m /V f and its associated parameters in terms 
of liquid filled levels, (B) Ratio of t^/tf and its 
associated parameters in terms of liquid filled levels. 

(A) Ratio of t B /t f and its associated parameters in terms 

\ 

of liquid filled levels, (B) Ratio of a m /a g and its 
associated parameters in terms of liquid filled levels. 
(A) Ratio of L m /h and its associated parameters in terms of 
liquid filled levels, (B) Ratio of V m /V fin and its 
associated parameters in terms of liquid filled levels. 
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(H) Cryogenic Liquid Hydrogen Reorientation 
Activated by High Frequency Impulsive 
Acceleration of Geyser Initiation 
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Abstract 


The requirement to settle or to position liquid fluid over the 
outlet end of spacecraft propellant tank prior to main enqine restart 
poses a microgravity fluid behavior problem. Resettlement or 
reorintation of liquid propellant can be accomplished by providing 
optimal acceleration to the spacecraft such that the propellant is 
reorientaed over the tank oulet without any vapor entrainment, any 
excessive geysering, or any other undesirable fluid motion for the 
space fluid management under microgravity environment. The purpose 
of present study is to investigate most efficient technique for 
propellant resettling through the minimization of propellant usage 
and weight penalties. Comparison between the constant reverse 
gravtiy acceleration and impulsive reverse gravity acceleration to be 
used for the activation of propellant resettlement, it shows that 
impulsive reverse gravity thrust is superior to constant reverse 
gravity thrust for liquid reorientation in a reduced gravity 
env ironment . 
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I . Introduction 


In spacecraft design, the requirements fora settled propellant 
are different for tank pressurization, engine restart, venting, or 
propellant transfer. Prepressurization requires that heat and mass 
transfer effects be minimized; otherwise, a process of chill down of 
tank, venting of noncondensing gases, etc. , may have to carry out for 
the cryogenic system. For engine restart, it is necessary to have the 
liquid settle with no bubbles near the tank outlet so that the initial 
flow of propellant will not carry vapor to the pump or engine. The 
slosh wave amplitude should be relatively low to keep the center of 
mass shifts within an acceptable range and wave motion low enough to 
avoid pressure collapse caused by interface agitation . For venting , 
it is probably necessary that virtually all bubbles be displaced from 
the bulk liquid so that a two-phase mixture is not vented. Propellant 
transfer requires that the liquid be completely settled with 
virtually no bubbles. Outflow of a liquid near the tank outlet can 
result in the prematrue ingestion of gas while a signigicant amount of 
liquid is still in the tank under microgravity environment. This 
phenomenon is termed "suction dip". Slosh wave motion must be 
minimal because the combination of "suction dip" and sloshing could 
cause gas pull-through to occur more readily in microgravity than if 

the surface were essentially quiescent. 

During the prepressurization of a cryogenic propellant in 
microgravity, significant heat and mass transfer will occur if the 
liquid interface is disturbed. Interface disturbances may result 
from (a) impingement of the gas on the liquid surface at a mass flow 
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rate sufficient to cause Kelvin -Helmholtz instability, (b) globule 
formation from breaking waves caused by wave motion over baffles or 
internal hardware, (c) globule and surface froth formation resulting 
from movement of bubbles through the liquid to the surface, and (d) 
surface froth formation because of gas impingement. 

Bubble and globule formation as a result of liquid impact with 
the aft end of the tank could lead to propellant loss for the 
spacecraft during venting. Globules could be entrained in the vented 
ullage gas or bubbles rising through the liquid and expanding because 
of the decreasing tank pressure could cause a spray of globules to be 
vented. Liquid level rise, vent liquid loss, fluid freezing, and 
vehicle dynamics are all affected by the microgravity levels. 

A key objective of the cryogenic fluid management of spacecraft 
propulsion system, such as a Space Transfer Vehicle 1 (STV) , is to 
develop the technology necessary for acquisition or positioning of 
liquid and vapor within a tank in reduced gravity to enable liquid 
outflow or vapor venting. Liquid acquisition techniques can be 
divided into two general categories: (1) Active liquid acquisition 
by the creation of a positive acceleration environment resulting from 
the propulsive thrust of small auxiliary engines, and (2) Passive 
liquid acquisition utilizing the liquid capillary forces provided by 
using solid baffles of liquid traps made of fine mesh screen material. 
In this series of study, active liquid acquisition is aimed for 
numerically simulating the resettlement of cryogenic liquid 
hydrogen. Liquid hydrogen, which, in general, poses more severe 
technical challenges than liquid oxygen, is used as the test bed 
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working fluid in this study. 

Recently Lesie 2 was able to measure and to numerically compute 
the bubble shapes at various ratios of centrifugal force to surface 
tension force in 2 , 4 and 6.3 cm deep cylinders in the microgravity 
environment. The results showed excellent agreement between model 
computation and measurements. Hung and Leslie 3 extended Leslie's 
work 2 to rotating free surfaces influenced by gravity with higher 
rotating speeds when the bubble intersects with both the top and 
bottom walls of the cylinder. Hung et al. 4 ' 5 further extended the 
work to include rotating speeds which resulted with bubbles 
intersecting and/or without intersecting the top, bottom and side 
walls of the cylinder. 

An analysis of time— dependent dynamical behavior of surface 
tension on partially— filled rotating fluids in both low gravity and 
microgravity environments was carried out by numerically solving the 
Navier-Stokes equations subjected to the initial and the boundary 
conditions 4 ' 6 . At the interface between the liquid and the gaseous 
fluids, both the kinematic surface boundary condition, and the 
interface stress conditions for components tangential and normal to 
the interface, were applied. The initial condition for the bubble 
profiles was adopted from steady-state formulations developed by Hung 
and Leslie 3 , and Hung et al. 6 for rotating cylinder tank; and by Hung 
et al. 5 ' 7 for the dewar-shaped container to be used in the Gravity 
Probe-B Spacecraft 8 . Some of the steady-state formulations of 
bubble shapes, in particular for bubbles intersecting at the top wall 
of the cylinder, were compared with the experiment carried out by 
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Leslie 2 in a free-fallinng aircraft (KC-135) . Comparisons of time- 
dependent results between numerical computations and experiments 
were unavailable. This was because the calibration of the recordings 
of time-dependent gravity variations in a KC-135 aircraft during the 
short time periods of microgravity environment is very difficult. 
There was also an unavailability of accelerometer data for measuring 
the actual levels of microgravity during the experiment. 

An efficient propellant settling technique should minimize 
propellant usage and weight penalties. This can be accomplished by 
providing optimal acceleration to the spacecraft such that the 
propellant is reoriented over the tank outlet without any vapor 
entrainment, any excessive geysering, or any other undesirable fluid 
motion . 

Production of geyser during the propellant reorientation is not 
a desirable motion for the space fluid management. It is because 
geyser is always accompanied by the vapor entrainment and globule 
formation. Geyser is observed at reverse gravity thrust greater than 
certain critical values of acceleration during the course of liquid 
reorientation. In other words, geyser will not be observed at very 
low reverse gravity level , and it will be detected when the reverse 
gravity level is greater than the certain critical value. In this 
series of study, numerical simulation of positive liquid acquisition 
is attemted by introducing reverse gravity acceleration, resulting 
from the propulsive thrust of small auxiliary engines which exceeds 
the critical value for geyser initiation. The reverse gravity 
acceleration is starting with a small value and increases gradually 
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till the initiation of geyser is detected in the computer simulation 
for the liquid reorientation of propellant tank with various liquid- 
filled levels. 

In this series of studies, time-dependent computations have been 
carried out to investigate the dynamical behavior of fluid 
reorientation or resettling of propellant prior to main engine firing 
for spacecraft restart at the net reverse gravity acceleration which 
is great enough to initiate geyser during the liquid reorientation. 
First paper of present study (Paper I ) 9 investigates the 
characteristics of fluid resettlement due to the reversal of lowest 
constant reverse gravity acceleration which is great enough to 
initiate geyser without any impulsive acceleration. The frequency 
of impulsive acceleration is generally termed "gravity jitters". 
Gravity j itters are produced by spacecraft attitude motion , machinery 
(turbine, pump, engine) vibrations, thruster firing, thruster 
shutdown, impulsive engine acceleration, etc . 10 Positioning of 
liquid propellant over the tank outlet can be carried out by using 
small auxiliary thrusters which provide a thrust parallel to the 
tank's major axis in the direction of flight. 

Computer simulation of flow field based on Paper I 9 during the 
course of fluid reorientation induced by constant reverse gravity 
acceleration show that six dimensionless parameters are abtained in 
the study. These parameters hold near constant values through the 
entire ranges of liquid filled levels during the course of fluid 
reorientation activated by the reverse gravity acceleration great 
enough to initiate geyser. As the denominators of these 
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dimensionless parameters are either predetermined from the geometry 
of liquid filled levels or can be deduced from the corresponding 
calculations associated with the geyser initiation gravity levels. 
One can predict the values of these flow parameters. These 
predictable parameters include maximum flow velocity V m , time for 
observing maximum flow velocity t m , time for reorienting liquid 
flowing down and reaching the bottom of propellant tank t R , scale 
length of maximum flow velocity L m , and scale flow acceleration 
associated with maximum velocity a m . 

Instead of applying constant reverse gravity acceleration as we 
described in Paper I 9 , this paper adopts impulsive reverse gravity 
acceleration with a high frequency of 10 Hz for the activation of fluid 
reorientation with liquid filled levels of 30, 50, 65, 70 and 80%. 

II. Mathematical Model 

The present study examines time-dependent fluid behaviors, in 
particular the dynamics of liquid hydrogen and hydrogen vapor 
reorientation induced by reverse gravity acceleration which is great 
enough to introduce geyser initiation. As in Paper I 9 , time- 
dependent axial symmetry mathematical formulation are adopted. 

Consider a closed circular cylinder of radius, a, with length, L, 
which is partially filled with a cryogenic liquid hydrogen of constant 
density p and kinematic viscosity v . Let us use cylindrical 
coordinates (r, 9, z) , with corresponding velocity components (u, v, 
w) . The gravitational acceleration, g, is along the z-axis. For the 
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case of axial symmetry, the 0-dependency vanishes. The governing 
equations are shown as follows: 


(A) Continuity Equation 

1 a 3w 

(ru) + — =0 (2-1) 

r 3r 3z 

(B) Momentum Equations 

Du v 2 1 3P ( u ' 

_ — = — — — + y 7 2 u (2-2) 

Dt r p dr l r 2 t 

Dv uv f 

— + — = y 7 2 v - 

Dt r l 

Dw 1 3P 

g + y 7 2 W (2-4) 

Dt P 3z 

where , 

D 3 9 3 

— = — + u — + w — (2-5) 

Dt 3t 3r 3z 


3 ) 3 2 

r — + (2-6) 

dr) 3z 2 


Let the profile of the interface between gaseous and liquid 
fluids be given by: 

t, (t, r, z) =0, or r = T(t,z) (2-7) 

initial condition of the profile of interface between 
gaseous and liquid fluids at t = t, is assigned explicitly, and is 
given by: 

t,(t = t 0 , r, z) =0, or r = n 0 (z) = v(t 0 ,z) (2-8) 

A set of boundary conditions has to be supplied for solving the 



v ' 
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equations. These initial interface profiles used in this study have 
been given explicitly through the steady state computations made by 
Hung and Leslie 3 and Hung et al 4 which were checked by the experiments 
carried overby Leslie 2 . These boundary conditions are as follows: 

(1) At the container wall, no-penetration and no-slip 
conditions assure that both the tangential and the normal 
components of the velocity along the solid walls will 
vanish. In the numerical calculation of bubble profiles 
for ethanol and air, a constant contact angle is present when 
the free surface of liquid ethanol intersects the container 
wall . 

(2) Along the interface between the liquid and gaseous fluids, 
the following two conditions apply: 

(a) Kinematic surface boundary condition: The liquid (or 

gaseous) surface moves with the liquid (or gas) which 
implies 


Dr, 


= 

0, or 

Dt 


dr) 

dr) dr) 

+ 

U + W 

at 

3r 9z 


on ^ =-o (t = tj , r, z) 

/ 

(b) Interface stress condition: At the interface, the stress must be 

continuous . These can be decomposed to the components normal 
and tangential to the interface. For the component tangential 


136 


to the interface between liquid and gaseous fluids 

r*n — (n*r*-n) = T*n — (n*r*n) (2-10) 

J liquid!. Jgas 

must hold. Here 

3 u i a u j 2 9 u k 9 u k 

T i i = U ( + + 5 i j ) + f 5 i j 

3 Xj 3 Xj 3 d x k 3x k 

is the viscous stress tensor; u, the viscous coefficient of the 
first kind? the viscous coefficient of the second kind; n, the 
unit vector normal to the interface and 5 { j , the Dirac delta 
function. For the component normal to the interface between the 
liquid and gaseous fluids, the expression becomes Laplace's 
formula which is 

Pq _ Pl - (n*T*n) gas + ( n * T * n ) 1 i q u i d 

( 2 - 11 ) 

Here P L denotes the liquid pressure at the interface; P G , the 
gaseous pressure at the interface; <r, the surface tension of the 
Interface; and i> , the tangent of the interface which is defined 
by: 

dz 

<f> = — on i)i = t|(t ir r r z) ( 2 - 12 ) 

dr 

Detailed description of computational algorithm applicable to 


<r d f Tt 

r dr ( 1 + t 2 ) 1 ' 2 
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microgravity fluid management are illustrated in Paper I 9 and our 
earlier studies . 4 ' 6 ' n ' 13 As we have indicated in Paper 1 9 , for the 
purpose of facilitating easy comparison between computational 
resuslts and experimental measurement, a model of 0.01 size prototype 
is adopted in the computer simulation. Model size is height L = 
4.23672 cm and diameter D = 4.2672 cm. If the spacecraft had been 
coasting for a long time, aligned with its direction of motion, the 
most signigicant force, drag, would be axial and with acceleration of 
10“ 4 g 0 along upward direction. The hydrogen vapor is, thus, 
originally positioned at the bottom of the tank. The requirement to 
settle or to position liquid fuel over the outlet end of the spacecraft 
propellant tank prior to main engine restart poses a microgravity 
fluid behavior problem. 11 Retromaneuvers of spacecraft, such as 
STV, require settling or reorientation of the propellant prior to main 
engine firing. 1,11 Cryogenic liquid propellant shall be positioned 
over the tank outlet by using auxiliary thrusters (or idle— mode 
thrusters from the main engine) which provide a thrust parallel to the 
tank's major axis in the direction of flight. Similar to Paper I 9 , a 
small value of reverse gravity acceleration (downward direction) is 
provided by the propulsive thrust of small auxiliary engine to 
initiate the reorientation of liquid propellant. This small value of 
reverse gravity acceleration of propulsive thrust increases 
gradually till reaching the critical value on which initiation of 
geyser is detected during the time period of fluid resettlement. We 
term this reverse gravity acceleration of propulsive thrust, which is 
capable to initiate geyser, as "geyser initiation gravity-level". 
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geyser ini.ti. 3 td. 0 n gravity level has been investigated through 
the method of trial and error for the various liquid-filled levels as a 
base to simulate impulsive reverse gravity acceleration with 
frequencies of 0.1, 1.0 and 10 Hz. As we have indicated in Paper I 9 , 
cryogenic liquid hydrogen at temperature of 2 OK is considered. 
Hydrogen density of 0.071 g/cm 3 ; surface tension coefficient at the 
interface between liquid hydrogen and hydrogen vapor of 1.9 dyne/ cm; 
hydrogen viscosity coefficient of 1.873 x 10 3 cm 2 /s; and contact 
angle of 0.5° are used in the computer simulation. 

III. Numerical Simulation of Liquid Hydrogen 
Reorientation With Geyser Initiation at High 
Frequency Impulsive Reverse Gravity 
Acceleration With Frequency of 10 Hz 

In this paper, among three categories of impulsive gravity 
acceleration with frequencies of 0.1, 1.0 and 10 Hz, reorientation of 
cryogenic liquid hydrogen activated by geyser initiation impulsive 
reverse gravity acceleration with a high frequency of 10 Hz produced 
by propulsive thrust will be investigated for various liquid filled 
levels of propellant tank. Paper I shows that these geyser 
initiation gravity levels are 5.5 x 10' 2 , 6.52 x 10 2 ,6.6xl0 2 ,6.7x 
10' 2 and 8.2 x 10" 2 g 0 for liquid filled levels of 30, 50, 65, 70, and 

80%, respectively. 

Table 1 shows some basic geometrys and characteristics of 
cryogenic liquid hydrogen resettlement activated by reverse gravity 
acceleration at geyser initiation gravity level. Average liquid 
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height h, and maximum liquid height h ra are shown in Figure 1. Average 
fall velocity V ^ ^ average free fall time t^ , and free fall 
velocity from maximum liquid height V fw are computed from the 
following equations: 


(2g i0 h) 1/2 

(3-1) 

(2h 1 

1/2 

(3-2) 

Wi<J 



(2g i( 

h ) 1/2 

j AA Tn / 

(3-3) 


where g io denotes geyser initiation reverse gravity acceleration. 

To show examples of the selected sequences of time evolution of 
fluid reorientation for cryogenic hydrogen. Figures 2, 3, 4, 5, and 6 
show time evolution of fluid reorientation activated by geyser 
initiation impulse reverse gravity acceleration with a frequency of 
10 Hz for liquid filled levels of 30, 50, 65, 70, and 80-s, 
respectively. Each figure contains four sub-figures. Subfigure 
(A) is initial profile of liquid-vapor interface at the moment of the 
starting of fluid reorientation at time t = 0; subfigure (B) , the flow 
profile during the course of fluid reorientation before the 
initiation of geysering motion; subfigure (C) , the flow profile with 
geysering motion; and subfigure (D) , the flow profile after the ending 
of geysering motion. 

Examples of selected sequences of time evolution of fluid 
reorientation illustrate following flow behaviors: (1) The liquid 
starts to flow in an annular sheet along the solid wall of tank and 
gradually pushes the vapor toward the central portion of the lower 
dome of tank as the net acceleration, reversing the direction of 
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gravity field, which is applied toward the downward direction of the 
tank's major axis, by using small auxiliary thrusters; (2) As the 
downward fluid annular sheet along the tank wall reaches the central 
bottom dome side of the tank, a geysering flow is observed; and (3) The 
vapor is thus pushed upward centrally into the liquid and the 
geysering disappears. 

Based on the computer simulation of flow field values of maximum 
flow velocity V TO , time for observing maximum flow velocity t m , and 
time for reorienting liquid flowing down and reaching the bottom of 
propellant tank t R , are obtained and illustrated in Table 2 for 
reverse gravity acceleration with impulsive frequency of 10 Hz. 
Scale length of maximum flow velocity L m , and scale flow acceleration 
associated with maximum velocity a m , can be computed from the 
following parameters: 


= v m t m 


a 


m 



(3-4) 

(3-5) 


Following dimensionless parameters are introduced: V m /V f/ t R /t f , 
t,/t f , a m /a g , L m /h and V m /V fm where a g stands geyser initiation 
acceleration (cm/s 2 ) for corresponding geyser initiation gravity 
level 9 i o * Impulsive reverse gravity acceleration, g { with 
frequency f Hz is defined as follows: 


9i 



sin 


\ 

2 TTf t 

/ 


(3-6) 


Figures 7 to 9 show the variations of dimensionless parameters in 
terms of liquid filled levels for impulsive acceleration with 
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Denominators of these six dimensionless 


frequency of 10 Hz. 
parameters are either predetermined from the geometry of liquid fill 
levels or can be deduced from the corresponding calculations 
associated with the geyser initiation gravity levels. 
Characteristics of these near constant values dimensionless 
parameters can provide a good understanding of the physics of 
microgravity fluid behaviors, in particular the active category of 
liquid acquisition or positioning, and also the design criteria of on- 
orbit spacecraft propulsion system at the critical value of reverse 
gravity acceleration of propulsive thrust which is capable to 
initiate geyser. 

An efficient propellant settling technique should minimize 
propellant usage and weight penalties. This can be accomplished by 
providing optimal acceleration to the spacecraft such that the 
propellant is reoriented over the tank outlet without any vapor 
entrainment, any excessive geysering, or any other undesirable fluid 
motion. In particular, it is important to study how well the 
impulsive acceleration can provide higher efficient propellant 
settling technique than the constant acceleration thrust technique. 
Also, what is the most optimal choice of impulsive frequency which can 
achieve the best fluid acquisition management being the goal of our 
research. 

Figure 7 (A) shows the ratio of maximum flow velocity to average 
free fall flow velocity V m /V f and its associated parameters of V m and 
V in terms of liquid filled levels for impulsive acceleration with 
frequency of 10 Hz. It shows that the ratio of V B /V f varies in the 
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range of 5.1 to 5.3 in the entire liquid filled levels while V m andV f 
vary from 63.0 to 89.3 cm/s (decreasing with increasing liquid filled 
levels) and from 12.5 to 17.0 cm/s (also decreasing with increasing 
liquid filled levels) , respectively. As V f can be predetermined from 
geyser initiation gravity level and average liquid height, shown in 
Equation (3-1) , one can make an approximate prediction of maximum flow 
velocity during the liquid reorientation for the various liquid 
filled levels. In comparison for impulsive frequency acceleration 
and constant reverse thrust acceleration in terms of the performance 
of capability to induce higher maximum flow velocity to accomplish 
flow reorientation, it shows that impulsive frequency of 10 Hz reverse 
thrust is better than that of constant thrust acceleration in which 
the maximum flow velocity induced are in the range of 62.0 to 73 cm/s in 

the entire liquid filled levels. 

Figure 7(B) shows the ratio of liquid reaching bottome time to 
average free fall time t R /tp and its associated parameters of t p and t^ 
in terms of liquid filled levels for impulsive acceleration with 
frequency of 10 Hz. It shows that the ratio of tp/t^ varies in the 
range of 1.1 to 1.2 in the entire liquid filled levels while t R and t f 
vary from 0.18 to 0.37 s (decreasing with increasing liquid filled 
levels) , respectively. As t f can be predetermined from geyser 
initiation gravity level and average liquid height, shown in Equation 
( 3 — 2 ) , one can predict the time for reorienting liquid fluid flowing 
down from the original position and reaching the bottom of propellant 
tank for the various liquid filled levels at the reverse gravity 
acceleration capable for the intiation of geyser. In comparison for 
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impulsive frequency of 10 Hz reverse thrust acceleration and constant 
thrust acceleration, in terms of the performance of capability to 
induce the flow which can reach the tank bottom with a minimum time, it 
shows that the performance for impulsive acceleration with 10 Hz is 
better than that of constant thrust acceleration in which t R is in the 
range of 0.20 to 0.40 s for the entire liquid filled levels. 

Figure 8 (A) shows the ratio of time for observing maximum flow 
velocity to average free fall time t m /t f and its associated parameters 
of t m andt f in terms of liquid filled levels. It shows that the ratio 
of t m /t f is in a value of 1.2 in the entire liquid filled levels while 
t m and t f vary from 0.19 to 0.39 s (decreasing with increasing liquid 
filled levels) and form 0.15 to 0.32 s (also decreasing with 
increasing liquid filled levels) , respectively. As we indicated in 
Figure 7(B) , t f can be predetermined, one can predict the time for 
observing maximum flow velocity for various liquid filled levels at 
the reverse gravity acceleration capable for the initiation of 
geyser. In comparison for impulsive frequency of 10 Hz reverse 
thrust acceleration and constant thrust accelerattion, in terms of 
performance of capability to imduce maximum flow velocity in a short 
period of time, it shows that the perf romance for impulsive 
acceleration with 10 Hz is better than that of constant thrust 
acceleration in which t m is in the range of 0.2 0 to 0.4 2 s for the 
entire liquid filled levels. 

Figure 8(B) shows the ratio of scale flow acceleration 
associated with maximum velocity to geyser intiation acceleration for 
corresponding geyser intiation gravity level a m /ag and its associated 
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parameters of a m and a g in terms of liquid filled levels for impulsive 
aceeleration with frequency of 10 Hz. It shows that the ratio of 
a m /a g vary from 4.2 to 4.3 in entire liquid filled levels while a m and 
a g vary from 228.9 to 331.6 cm/s 2 (increasing with increasing liquid 
filled levels) and from 53.9 to 80.4 cm/s 2 (also increasing with 
increasing liquid filled levels), respectively. As a g can be 
predetermined from geyser initiation gravity level, one can make an 
approximate prediction of scale flow acceleration associated with 
maximum velocity, which is defined in Equation (3-5) , at the reverse 
gravity acceleration capable for the initiation of geyser. In 
comparison for impulsive frequency of 10 Hz reverse thrust 
acceleration and constant thrust acceleration, in terras of the 
performance of capability to induce higher scale flow acceleration 
associated with maximum flow velocity, it shows that the performance 
for impulsive acceleration with 10 Hz is better than that of constant 
thrust acceleration in which a m is in the range of 175 to 310 cm/s 2 for 
the entire liquid filled levels. 

Figure (9A) shows the ratio of scale length of maximum flow 
velocity to average liquid height L m /h and its associated parameters 
of L m and h in terms of liquid levels for impulsive acceleration with 
frequency of 10 Hz. It shows that the ratio of L TO /h varies in the 
range of 12 . 8 to 13 . 0 in the entire liquid filled levels while L m and h 
vary from 11.9 to 34.8 cm (decreasing with increasing liquid filled 
levels) and from 0.93 to 2.67 cm (also decreasing with increasing 
liquid filled levels) , respectively. As h can be predetermined from 
the geometry of liquid filled levels, one can make an approximate 
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prediction of scale length of maximum flow velocity , which is defined 
in Equation (3-4) , at the reverse gravity acceleration capable for the 
initiation of geyser. In comparison for impulsive frequency of 10 Hz 
rsv6rs6 thrust acceleration and constant thrust acceleration, in 
terms of the performance of capability to produce larger scale length 
of maximum flow velocity, it shows that the scale length produced by 
impulsive acceleration with 10 Hz is longer than that of constant 
thrust acceleration in which L m is in the range of 12.6 to 3 0.9 cm for 
the entire liquid filled levels. 

Figure 9 (B) shows the ratio of maximum flow velocity to free fall 
velocity from maximum liquid height V m /V fTn and its associated 
parameters of V m and V f m in terms of liquid filled levels for impulsive 
acceleration with frequency of 10 Hz. It shows that the ratio of 
V /v^ varies in the range of 3.9 to 4.6 in the entire liquid filled 
levels while V m and V fTn vary from 63.0 to 89.3 cm/s (decreasing with 
increasing liquid filled levels) and from 16.4 to 19.2 cm/s (also 
decreasing with increasing liquid filled levels) , respectively. As 
V fin can be predetermined from geyser initiation gravity level and 
maximum liquid height, shown in Equation (3—3) , one can make an 
approximate prediction of maximum flow velocity at the reverse 
gravity acceleration capable for the initiation of geyser. In 
comparison for impulsive frequency of 10 Hz reverse thrust 
acceleration and constant thrust acceleration, in terms of producing 
higher maximum flow velocity, it shows that the maximum flow velocity 
produced by impulsive acceleration with 10 Hz produce higher maximum 
flow velocity than that of the constant thrust acceleration in which 
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V m is in the range of 62.0 to 73.6 cm/s for the entire liquid filled 
levels. 

Similar to Paper I 9 for constant reverse gravity acceleration, 
six dimensionless parameters presented in this study show that the 
parameters hold near constant values through the entire ranges of 
liquid filled levels during the course of reorientation of liquid 
hydrogen activated by the reverse gravity acceleration with impulsive 
thrust frequency of 10 Hz, which is great enough to initiate geysering 
flow. As the purpose of present study is to investigate an efficient 
propellant resettling technique which is able to minimize propellant 
usage and weight penalties, comparison of flow parameters induced by 
the impulsive acceleration with frequency of 10 Hz and constant thrust 
acceleration toward these goals have been made. It shows that the 
operation of auxiliary engine with impulsive thrust acceleration of 
10 Hz is better than that of the constant thrust acceleration. 

IV. Discussion and Conclusions 

The requirement to settle or to position liquid fluid over the 
outlet end of the spacecraft propellant tank prior to main engine 
restart poses a microgravity fluid behavior problem. Retromaneuvers 
of spacecraft require settling or reorientation of the propellant 
prior to main engine firing. Cryogenic liquid propellant is 
positioned over the tank outlet by using small auxiliary thrusters (or 
idle-mode thrusters from the main engine) which provide a thrust 
parallel to the tank's major axis in the direction of flight. 

The results of the study of fluid reorientation have to be 
evaluated in terms of how well they can be managed efficiently. An 
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efficient propellant settling technique should minimize propellant 
usage and weight penalties through the operation of small thrusters 
(or idle-mode thrusters from the main engine) . This can be 
accomplished by providing optimal acceleration to the spacecraft such 
that the propellant is reoriented over the tank outlet without any 
vapor entrainment, any excessive geysering, or any other undesirable 
fluid motion. 

Production of geyser during the propellant reorientation is not 
a desirable motion for the space fluid management under microgravity 
environment. It is because geyser is always accompanied by the vapor 
entrainment and globule formation. Geyser is observed at reverse 
gravity acceleration greater than certain critical values of 
acceleration during the course of liquid reorientation. In this 
paper, numerical simulation of positive liquid acquisition is 
attempted by introducing reverse gravity acceleration, resulting 
from the propulsive thrust with high impulsive frquency of 10 Hz 
auxiliary engine, which exceeds critical value for the initiation of 
geyser. 

Evaluation of performance is based on how efficient the 
impulsive reverse gravity with 10 Hz frequency in comparison with 
constant thrust acceleration can activate following flow parameters 
at the same background thrust accelerations: (A) a higher maximum 
flow velocity, (B) a shorter time period for flow to reach maximum 
velocity, (C) a shorter time period for flow to reach tank bottom (tank 
outlet) for fluid resettling, (D) a larger length scale of maximum 
flow velocity, and (E) a higher flow acceleration associated with 
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maximum flow velocity. Comparison between the results of present 
study for impulsive thrust with 10 Hz frequency and that of constant 
thrust, shown in Paper I 9 , it shows that impulsive thrust is superior 
than the constant thrust in terms of efficient operation of fluid 
reorientation . 

Based on the computer simulation of flow fields during the course 
of fluid reorientation, six dimensionless parameters are presented 
both in Paper I 9 and this study. It is shown that these parameters 
hold near constant values through the entire ranges of liquid filled 
levels during the course of fluid reorientation activated by the 
reverse gravity acceleration in both constant and impulsive thrusts 
great enough to initiate geyser. As the denominators of these 
dimensionless parameters are either predetermined from the geometry 
of liquid filled levels, as shown in Table 1, or can be deduced from the 
corresponding calculations associated with the geyser initiation 
gravity levels, one can predict the values of these flow parameters. 
Present study can greatly enhance our understanding in the behaviors 
of cryogenic fluid resettlement under reduced gravity environment. 
This is particularly important for liquid acquisition technique to be 
used in on-orbit spacecraft design. 

Any fluid capable of motion relative to the spacecraft will be 
subject to an acceleration relative to the mass center of the 
spacecraft that arises from the gravity gradient of the Earth 14 ' 1S . 
In addition to the Earth's gravitational force, the interaction 
between the particle mass of fluids and the spacecraft mass due to 
gravity gradient acceleration 1 4 have also been taken into 


149 


consideration in this microgravity fluid management study. 

To conclude, we have demonstrated that, the computer algorithm 
presented, can be used to simulate fluid behavior in a microgravity 
environment, in particular the development of technology necessary 
for acquisition or positioning of liquid and vapor within a tank to 
enable liquid outflow or vapor venting through active liquid 
acquisition by the creation of a positive acceleration environment 
resulting from propulsive thrust. Better understanding of the full 
pictures of flow fields in both constant and impulsive thrusts, during 
the course of fluid reorientation can provide the proper design 
techniques for handling and managing the cryogenic liquid propellants 
to be used in on-orbit spacecraft propulsion. It is important to 
emphasize that impulsive reverse gravity thrust is superior to 
constant reverse gravity thrust for the activation of liquid 
necessary for the resettlement of liquid in a reduced gravity 
environement . 
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TABLE 1 SOME BASIC GEOMETRYS AND CHARACTERISTICS 
OF CRYOGENIC LIQUID HYDROGEN REORIENTATION 
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Table 2 CHARACTERISTICS OF CRYOGENIC HYDROGEN REORIENTATION 
(Frequency of Impulsive Acceleration * 10 Hz) 
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Figure Captions 


Figure 1. 


Figure 


Figure 


Figure 


Figure 


(A) Distribution of grid points in the radial— axial plane 
of cylindrical coordinate for propellant tank, and (B) 
Model size propellant tank adopted for numerical 
simulation with geometrical description. 

. Selected sequences of time evolution of fluid 

reorientation with liquid filled level of 30% for 

impulsive acceleration with frequency of 10 Hz, (A) 

initial profile, (B) flow profile before the initiation of 
geyser, (C) flow prof ile with geyser, and (D) flow profile 
after the ending of geyser. 

3 . Selected sequences of time evolution of fluid 

reorientation with liquid filled level of 50% for 

impulsive acceleration with frequency of 10 Hz, (A) 

initial profile, (B) flow profile before the initiation of 
geyser, (C) flow prof ile with geyser , and (D) flow prof ile 
after the ending of geyser. 

4. Selected sequences of time evolution of fluid 

reorientation with liquid filled level of 65^ for 

impulsive acceleration with frequency of 10 Hz, (A) 

initial profile, (B) flow profile before the initiation of 
geyser, (C) flow prof ile with geyser , and (D) flow prof ile 
after the ending of geyser. 

5. Selected sequences of time evolution of fluid 

reorientation with liquid filled level of 70% for 

impulsive acceleration with frequency of 10 Hz, (A) 
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Figure 


Figure 7 . 


Figure 8 . 


Figure 9 . 


initial profile, (B) flow prof ile before the initiation of 
geyser, (C) flow profile with geyser, and (D) flow profile 
after the ending of geyser. 

. Selected sequences of time evolution of fluid 
reorientation with liquid filled levels of 80% for 
impulsive acceleration with frequency of 10 Hz, (A) 
initial profile, (B) flow profile before the initiation of 
geyser, (C) flow prof ile with geyser , and (D) flow profile 
after the ending of geyser. 

(A) Ratio of V m /V f and its associated parameters in terms 
of liquid filled levels, (B) Ratio of t R /t f and its 
associated parameters in terms of liquid filled levels. 
(A) Ratio of t^/tf and its associated parameters in terms 
of liquid filled levels, (B) Ratio of a m /a^ and its 
associated parameters in terms of liquid filled levels. 
(A) Ratio of L m /h and its associated parameters in terms of 
liquid filled levels, (B) Ratio of V m /V fm and its 
associated parameters in terms of liquid filled levels. 


156 






Liquid Hydrogen and Vapor 
Liquid Fllled-307. 1- 0-“ U3 

g — 5. 50 H l O'* g, 

r- i.oomiq'hz 


Liquid Hydrogen and Vapor 
Liquid Fllled-30X t“ 3.60*10 s 
g— 5.50 h10' ? 9i 
f- l .00* 10* Hz 


7Z 

t 5 

/ \ 

A 

a 

t 

' 1 


t— _ 

t 


\ 

X o 

1 


\ 

CD cm ” 

i 

1 

1 

' LJ 

1 1 

l / 

/ 

31 - 


■/ 


- 1.0 0-0 i -0 

RADIUS (CM) 


-l.o o.o i.o 
RADIUS (CM) 


Liquid Hydrogen and Vapor 
Liquid Fllled-30X 5.20*10 ' s 
g”~5.S0*l0 


r- i.oo*iq Hz 



Liquid Hydrogen and Vapor 
Llqutd Fulled“30X L“ 8.80*10 s 

g— 5.S0H lO'^g, 
f- 1.00*l0 l Hz 


RADIUS (CM) 


RAO I US (CM) 



HEIGHT I CN) HEIGHT (CHI 

2.0 3.0 4.0 r\ p.O 1.0 2*0 3.0 


Liquid Hydrogen and Vapor 

LLquLd Flll9d“50X 0.00s 
g—~6. 52*10 *g, 
f- 1.00*10 l Hz 


Liquid Hydrogen and Vapor 

liquid Fllled-507. 2.40*10 s 
g— 6.52*10 *g, 
f- 1.00*lQ l Hz 




-l.a 0.0 i.o 
RflOIUS (CM) 


- 2.0 - 1.0 


RflOIUS (CM) 


Liquid Hydrogen and Vapor 
Liquid FLlled-50Z t« 3.43*10'' s 
g- -6. 52*10 J g # 

r- 1 .00*1 0 l Hz 


D. Liquid Hydrogen and Vapor 
Liquid Fllled“50X t* 5.23*10 s 
g — 6.52*10 J g, 
f- L.00*l0 l Hz 


//Mir 

/ t i \ \ \ 

y / / \ \ \ 


RflOIUS (CM) 


RflOIUS (CM) 



HEIGHT (CM) HEIGHT (CM) 

2.0 3.0 4.0 f-v 0.0 i.o 2.0 3.0 


Liquid Hydrogen and Vapor 

LLquLd PLLled“65X t- 0.00s 
g—~6. 60*10 *g # 

f- l.OO*lO l Hz 



B. Llq uld Hydrogen and Vapor 
LLquLd FLLled-657 t- 2.40*10's 
g— ~ 6.60x10 2 g # 
f- 1 .00x10* Hz 



RFID I US (CM) 


RflOIUS ( CM 3 


Liquid Hydrogen and Vapor 
LLquLd FLlled-657 t- 3.37x10'' s 
g- -6. 60 *10 1 g t 
f- 1 .00x1 0* Hz 


D. Liquid Hydrogen and Vapor 
■LLquLd FLllad-657. t* 7.57*l0's 

g- - 6.60*10 J g^ 
f- l.OOxlO'Hz 





RflOIUS (CM) 


RflOIUS (CM) 



HEIGHT (CM) HEIGHT (CM) 

2.0 3.0 V° 0.0 1.0 2.0 3.0 


Liquid Hydrogen and Vapor 
LLcjuLd FUAed-707 t- 0.00s 
g--7. 00*10 J g t 


f- 1.00*10 Hz 


Liquid Hydroger and Vapor 
Llqutd Fllled-707 t- 1.80x10'' s 
g— 7.00*10 J g t 
f- I.OOhIO'Hz 




o.o i. 
RADIUS (CM) 


RAO I US (CM) 


Liquid Hydrogen ar d Vapor 
LLquLd FLlled-707. t- 3.33x10'' s 
g--7.G0*10’ J g, 
f- 1 .00* 1 0* Hz 


p. Liquid Hydrogen and Vapor 
LLquLd FLlled“70X L* 8-73*10 s 

g--7.00*10 J g, 

f- 1.00*10 ‘Hz 



RADIUS (CM) 


RADIUS (CM) 





Fig 7 (A, B) 











Fig 9 


MflXinUII LIQUID HEIGHT,7( m (cm) VELOClTY,L„-V.,U (cm) 





References 


( 1) NASA Office of Aeronautics and Space Technology , Technology for 

Future NASA Missions: Civil Space Technology Initiative and 

Pathfinder . NASA CP-3016, 1988. 

(2) Hung, R. J. , and Leslie F. W. , Bubble Shapes in a Liquid-Filled 

Rotating Container Under Low Gravity, Journal Spacecraft and 
Rockets , 25 . 70-74, 1988. 

(3) Hung, R. J. , Tsao, Y. D. , Hong, B. B. , and Leslie, F. W. , Time 
Dependent Dynamical Behavior of Surface Tension on Rotating 
Fluids Under Microgravity Environment, Advances in Space 
Research . 8., No. 12, 205-213, 1988. 

(4) Hung, R. J. , Tsao, Y. D. , Hong, B. B. , and Leslie, F. W. , 
Axisymmetric Bubble Profiles in a Slowly Rotating Helium Dewar 
Under Low and Microgravity Environments, Acta Astronaut ica . , . 
19 . 411-426, 1989. 

(5) Hung, R. J. , Tsao, Y. D. , Hong, B. B. , and Leslie, F. W. , Bubble 
Behaviors in a Slowly Rotating Helium Dewar in Gravity Probe-B 
Spacecraft Experiment , J . Space craft and Rockets j_ 2 6 U 167 — 17 2 , 
1989 . 

(6) Hung, R. J., Tsao, Y. D. , Hong, B. B. , and Leslie, F. W. , 
Dynamical Behavior of Surface Tension on Rotating Fluids in Low 
and Microgravity Enviroments, Interna tional J ournal for 
Mi croaravitv Research and Applications , 11, 81-95, 1989. 

(7) Hung, R. J., Lee, C. C. , and Tsao, Y. D. , Fluid Behavior in 
Microgravity Environment, Proceed . National Sc i .. Counc i l ( A ) , 


157 


14 . 21-34, 1990. 

(8) Hung, R. J., Tsao, Y. D. , and Leslie, F. W. , Effect of G-Jitters 
on the Stability of Rotating Bubble Under Microgravity 
Environment, Acta Astronautica . 21 , 309-321, 1990. 

(9) Hung, R. J. , Tsao, Y. D. , Hong, B. B. , and Leslie, F. W. , 

Dynamics of Surface Tension in Microgravity Environment, 

Progress in Aeronautics and Astronautics , 127 , 124-150 , 1990. 

(10) Hung, R. J., Lee, C. C. , and Shyu, K. L. , Reorientation of 

Rotating Fluid in Microgravity Environment With and Without 
GT-avity .Titters. J. Spacecraft and Rockets . 27, in press, 1990. 

(11) Hung, R. J. , Super fluid and Normal Fluid Helium II in a Rotating 
Tank Under Low and Microgravity Environments, Proceed. 
National Sci. Council (A) . 14, 289-297, 1990. 

(12) Hung, R. J. , Lee, C. C. , and Leslie, F. W. , Response of Gravity 
Level Fluctuations on the Gravity Probe-B Spacecraft 
Propellant System, J. Propulsion and Power , 6, in press, 1990. 

(13) Hung, R. J-, Lee, C. C. , and Leslie, F. W. , Gravity Jitter 

Response Slosh Wave Excitation on the fluid in a Rotating Dewar, 
Advances in Space Research . 11, in press, 1990. 

(14) Hung, R. J., and Shyu, K. L. , Cryogenic Hydrogen Reorientation 
and Geyser Initiation at Various Liquid-Filled Levels in 
Microgravity, Advances in Space Research . 11 , in press, 1990. 

(15) Hung, R. J. , Lee, C. C. , and Leslie, F. W. , Gravity-Jitter 

Effected Slosh Waves on the Stability of Rotating Bubble Under 
Microgravity Environment, Advances in Space Research , li, in 
press, 1990. 


158 



